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Block  19sABSTRACT 

The  height  and  velocity  of  visible  clouds  ind . percent  cover  at  several 
altitudes  over  a  portion  (£8  degree  cone)  of  the  sky  can  been 
determined  from  the  ground  by  means  of  triangulat ion  with  an  instrument 
consisting  of  a  pair  of  CIO  cameras  and  a  computer.  A  system  has  been 
built  which  is  stand-alone,  automatic  and  produces  reports  every  ten 
minutes.  Results  are  repeatable  and  accuracy  as  determined  by  indirect 
calibration  is  within  ten  percent.  The  system  can  operate  at  night  and 
in  light  snow  or  rain.  The  cameras  can  distinguish  features  in  a  cloud 
cover  that  is  featureless  to  the  naked  fye.  Cost  of  the  instrument 
should  not  exceed  *10,000. 
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FINAL  REPORT 

Th*  Whole  Sky  Sensor  Project — Phase  On* 
ABSTRACT 


he  height  and  velocity  of  visible  clouds  and  percent  cover  at  several 
altitudes  over  a  portion  <28  degree  cone)  of  th*  sky  can  been 
determined  from  the  ground  by  means  of  triangulation  with  an  instrument 
consisting  of  a  pair  of  CIO  cameras  and  a  computer.  A  system  has  been 
built  which  is  stand-alone,  automatic  and  produces  reports  every  ten 
minutes.  Results  are  repeatable  and  accuracy  as  determined  by  indirect 
calibration  is  within  ten  percent.  The  system  can  operate  at  night  and 
in  light  snow  or  rain.  Th*  cameras  can  distinguish  features  in  a  cloud 
cover  that  is  featureless  to  the  naked  eye.  Cost  of  th*  instrument 
should  not  exceed  *10,808. 

A  simple  algorithm  for  th*  corresponding  of  points  in  3-space  has  been 
found.  There  is  some  indication  that  types  of  clouds  can  be  identified 
from  an  examination  of  th*  population  statistics  of  th*  images. 

The  system  should  include  shelters  for  th*  instruments  and  an  auto-iris 
attachment.  More  experience  with  different  kinds  of  clouds  and 


extensive  field  testing  are  required, 
viewing,  non-hori*ontal  cloud  forms  and 
low  altitude  images  should  be  developed. 
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1.  Objectives  of  the  study 

The  purpose  of  this  project  mss  to  construct  an  instrument  consisting 
of  two  electronic  cameras  and  a  computer,  to  test  its  utility  and 
effectiveness  in  determining  the  height  and  velocity  of  clouds  in  the 
sky  and  per  cent  cover,  and  to  test,  modify  and  as  required  add  to  a 
collection  of  computer  algorithms  which  dealt  with  correspondence, 
pattern  recognition  and  interpretat ion  in  the  simulation  mode. 


2.  Description  of  the  instrument 


2. 1  Hardware! 

The  instrument  consists  of  two  CID  cameras  mounted  963.52  feet  apart  at 
the  same  elevation,  240.0  feet  PSL,  controlled  by  and  transmitting  data 
via  coaxial  cable  to  a  computer  which  analyzes  the  data  and  produces 
reports  every  ten  minutes  of  the  height  of  visible  layer*  of  clouds, 
their  velocity  ano  direction  and  per  cent  cover. 

Each  camera  is  mounted  in  a  gimballed  frame  and  leveled  to  point 
straight  up.  It  is  also  aligned  so  that  its  x-axis  passes  through  the 
center  of  the  other  camera.  The  measured  line  of  the  x-axis  lies  N 
62o25’  04 "W  true.  The  survey  was  don*  by  the  licensed  firm  of 
Chamberlin,  Wolford  and  Chin  of  Leesburg,  VP.  P  map  and  details  of  the 
installation  are  given  as  Figure  1. 

The  cameras  are  General  Electric  Model  TN2505  surveillance  cameras  each 
employing  a  360  by  240  CID  chip  and  several  features  which  permit  time 
exposures  and  auto-iris  control. 

Three  75-ohm  coaxial  lines  and  60  htz  110  volt  power  are  connected  to 
each  camera.  The  coaxial  lines  carry  the  analog  signal  (equivalent  to 
a  cable  TV  signal)  to  the  computer,  a  clock  signal  for  synchronizing 
the  cameras  with  the  computer,  and  control  line  for  operating  the 
" inject-inhibit"  (time  exposure)  feature.  P  conventional  converter  at 
each  camera  converts  110  v  house  current  to  12  v.  dc  for  the  operation 
of  the  cameras.  The  auto  iris  feature  is  not  implemented. 

The  lenses  are  conventional  IS  mm  fl6-fl.7  Fujinon  TV  lenses.  No 
optical  (light)  filters  have  been  tested,  although  the  spectral 
response  of  the  CID  chips  extends  well  into  the  infra-red  so  that 
additional  information  is  available  by  . selection  of  that  portion  of  the 
spectrum. 

In  the  computer,  two  boards,  one  for  each  camera,  are  installed  to 
interface  the  cameras  to  the  computer.  These  boards  or 
"fram^-grabbers"  are  manufactured  by  Epix,  Inc.  Each  contains  a  flash 
analog-to-digital  converter  which  converts  the  analog  TV  signa.  from 
the  camera  into  a  series  of  86400  8-bit  numbers  <0. .255)  per  frame  and 
stores  these  in  a  1  megabyte  ROM  30  times  per  second.  The  boards  also 
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output  the  received  signal  to  a  standard  TV  monitor  so  that  what  the 
camera  sees  can  be  seen  by  the  operator.  After  the  frame-grabbers 
collect  two  complete  images  from  each  camera  2  seconds  apart,  the 
computer  transfers  the  data  from  the  frame-grabber  ram  to  computer  ram 
for  processing. 


The  computer  is  a  PC’s  Limited  28&-12  AT  personal  computer  operating  at 
12  Megahertz,  about  the  fastest  inexpensive  computer  available.  It 
features  a  20-megabyte  hard  disk,  an  80288/80287  chip  set  and  1 
megabyte  of  RAM.  A  number  of  other  features  required  for  program 
development  are  included,  but  are  not  necessary  for  automated 
instrument  operation  itself. 

Specifications  of  the  cameras,  frame  grabbers  and  computer  are  given  in 
Annex  1 . 


2.2  System  software: 

System  software  consists  of  a  standard  MS-DOS  package,  Epix’s  SVPCIP 
system  which  controls  and  permits  access  to  the  frame  grabbers,  the 
Run-C  Professional  C  interpreter  (RCP),  used  to  edit  and  test  program 
modules  and  function  in  interpretive  mode,  and  the  Lattice  C  compiler 
which  is  used  to  compile  into  machine  language  programs  written  and 
debugged  in  RCP.  All  of  these  have  been  installed  in  the  computer  and 
used  as  required  during  the  development  process.  In  a  field 
instrument,  portions  of  these  will  be  transferred  from  their 
environments  and  compiled  into  a  single  machine  language  program 
eliminating  the  need  for  keyboard  management,  the  current  mode  of 
operation.  Speci f icat ions  of  these  systems  are  given  in  Annex  2. 

2U  3  Application  software: 

The  programs  written  by  the  Principal  Investigator  for  this  instrument 
are  in  three  groups.  The  first  group  of  routines  controls  the 
acquisition  of  cloud  Images,  turning  the  cameras  on,  taking  single 
frame  pictures,  and  storing  part  of  each  of  the  four  pictures  into  a 
single  file  on  hard  disk  so  that  it  can  be  shown  on  the  monitor  if 
required. 

The  second  group  calculates  the  altitudes  of  the  cloud  layers  from  the 
images,  then  their  velocities,  then  converts  these  numbers  into  a 
report  which  is  displayed  on  the  screen,  and  can  be  printed  or 
transmitted  by  modem. 


The  third  group  is  a  collection  of  service  routines  which  provide  for 
£*  display  of  the  various  images,  graphing  of  results,  lens  compensation, 

^  and  smoothing  the  signals  by  means  of  a  fast-Fourier  algorithm,  and 

.BAT  procedures  for  data  management.  This  last  group  is  not  part  of 
}u  the  operating  system,  but  was  written  to  assist  in  its  development, 

k-  The  programs  are  listed  and  described  in  detail  in  Annex  3. 


I 


2.4  Operation: 

At  start-up,  the  image  is  displayed,  irises  and  inject-inhibit  set  for 
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the  proper  exposure,  and  automatic  operation  started.  In  normal 
operation,  the  instrument  runs  all  the  time.  A  clock  in  the  computer 
is  queried  continuously  for  the  end  of  the  inter-observation  interval 
(usually  10  minutes). 

When  this  time  is  reached,  the  image  is  sampled,  the  end  of  the  next 
such  interval  is  calculated  and  stored,  and  the  end  of  the  inter-sample 
interval  (usually  2  seconds)  is  calculated  and  stored.  Then  a 
single-frame  picture  is  taken  by  each  camera  simultaneously,  converted 
by  that  camera’ s  ADC  and  stored  in  grabber  ram  (QRAM). 

At  the  end  of  the  inter-sample  interval,  another  pair  of  pictures  is 
taken  and  stored  at  a  different  location.  Then  QRAM  is  copied  to 
computer  ram  and  the  two  pairs  of  pictures  analyzed,  producing  the 
layers  and  the  velocities  of  the  layers  and  an  estimate  of  cloud  cover. 
When  analysis  is  complete,  the  report  is  prepared  and  transmitted  or 
displayed  and  the  program  waits  until  the  end  of  the  inter-observat ion 
interval,  at  which  time  the  process  repeats. 

In  this  Phase  I  development  the  procedure  has  been  to  set  manually  the 
irises  on  *he  basis  of  the  pictures  being  displayed  continuously  on  the 
monitor.  When  the  exposure  is  correct,  the  automatic  operation  begins. 
With  the  exception  of  auto-iris,  which  is  not  implemented  in  this 
system,  the  automatic  system  is  identical  to  the  .  developmental  system 
at  this  time. 


2. S  Proceduresi 

The  procedures  divide  into  system  control,  data  management,  overlaying 
algorithms  and  report  preparation. 

System  control  consists  of  querying  the  realtime  clock  and  executing 
the  proper  routines  at  the  proper  time  as  described  above,  and  includes 
start  up. 

Data  management  is  more  complicated  for  two  reasons.  First,  the 
largest  data  object  that  can  be  handled  easily  by  the  system  software 
is  limited  to  65536  bytes,  while  the  image  produced  by  the  cameras  is 
86400  bytes.  Hence  the  image  must  be  dealt  with  in  chunks.  Similarly, 
the  "window"  from  grabber  RAM  into  computer  RAM  is  also  65536  bytes, 
thus  each  image  must  be  partitioned,  copied  21600  bytes  (the  most 
convenient  chunk)  at  a  time  into  computer  RAM,  compressed,  transformed 
to  integers,  and  filed  as  a  record  for  monitor  display  if  desired. 

Secondly,  because  of  optical  distortion  caused  by  the  fact  that  the 
lens  of  the  camera  is  rot  a  pinhole,  and  because  the  corners  of  the 
image  do  not  receive  the  same  light  intensity  as  the  center,  only  the 
center  300  columns  of  the  center  120  rows  are  usable  without 
substantial  compensat ion.  This  subset  is  extracted  from  the  raw  data 
and  refiled  as  arrays  of  integers  for  later  use.  Thus  the  data  goes 
through  a  number  of  steps  between  each  analytic  sequence,  and  a  number 
of  functions  are  devoted  entirely  to  this  activity. 
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the  proper  eKposure,  and  automatic  operation  started.  In  normal 
operation,  the  instrument  runs  all  the  time.  A  clock  in  the  computer 
is  queried  continuously  for  the  end  of  the  inter-observat ion  interval 
(usually  10  minutes). 

When  this  time  is  reached,  the  image  is  sampled,  the  end  of  the  next 
such  interval  is  calculated  and  stored,  and  the  end  of  the  inter-sample 
interval  (usually  P.  seconds)  ,  is  calculated  and  stored.  Then  a 
single-frame  picture  is  taken  by  each  camera  simultaneously,  converted 
by  that  camera’s  ADC  and  stored  in  grabber  ram  (GRAM). 

At  the  end  of  the  inter-sample  interval,  another  pair  of  pictures  is 
taken  and  stored  at  a  different  location.  Then  GRAM  is  copied  to 
computer  ram  and  the  two  pairs  of  pictures  analyzed,  producing  the 
layers  and  the  velocities  of  the  layers  and  an  estimato  of  cloud  cover. 
When  analysis  is  complete,  the  report  is  prepared  and  transmitted  or 
displayed  and  the  program  waits  until  the  end  of  the  inter-observat ion 
interval,  at  which  time  the  process  repeats. 

In  this  Phase  I  development  the  procedure  has  been  to  set  manually  the 
irises  on  the  basis  of  the  pictures  being,  displayed  cont inuous' y  on  the 
monitor.  When  the  exposure  is  correct,  the  automatic  operation  begins. 
With  the  exception  of  auto-iris,  which  is  not  implemented  in  this 
system,  the  automatic  system  is  identical  to  the  developmental  system 
at  this  time. 


2.5  Proceduresi 

The  procedures  divide  into  system  control,  data  management,  overlaying 
algorithms  and  report  preparation. 

System  control  consists  of  querying  the  realtime  clock  and  executing 
the  proper  routines  at  the  proper  time  as  described  above,  and  includes 
start  up. 

Data  management  is  more  complicated  for  two  reasons.  First,  the 
largest  data  object  that  can  be  handled  easily  by  the  system  software 
is  limited  to  65536  bytes,  while  the  image  produced  by  the  cameras  is 
66400  bytes.  Hence  the  image  must  be  dealt  with  in  chunks.  Similarly, 
the  "window"  from  grabber  RAM  into  computer  RAM  is  also  65536  bytes, 
thus  each  image  must  be  partitioned,  copied  21600  bytes  (the  most 
convenient  chunk)  at  a  time  into  computer  RAM,  compressed,  transformed 
to  integers,  and  filed  as  a  record  for  monitor  display  if  desired. 

Secondly,  because  of  optical  distortion  caused  by  the  fact  that  the 
lens  of  the  camera  is  not  a  pinhole,  and  because  the  corners  of  the 
i mage  do  not  receive  the  same  light  intensity  as  the  center,  only  the 
center  300  columns  of  the  center  120  rows  are  usable  without 
substantial  compensat ion.  This  subset  is  extracted  from  the  raw  data 
and  ref i led  as  arrays  of  integers  for  later  use.  Thus  the  data  goes 
through  a  number  of  steps  between  each  analytic  sequence,  and  a  number 
of  functions  are  devoted  entirely  to  this  activity. 
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time,  and  the  values  in  the  cells  of  one  subtracted  from  the  matching 
cells  of  the  other.  The  difference  for  each  cell  is  squared  and  summed 
and  the  total  divided  by  the  number  of  cells  which  were  differenced  for 
this  offset  (ofs).  In  computerese 

for  each  offset (ofs) |  efs»0.  .149 

for  each  row ( i ) |  i *0. . 59 

for  each  column(j);  j*0. .149 

del«zimg0(0, i,  )>-zimg0(l, i,  j+ofs) f 
sumsq*sumsq+del*del 
next  column 
next  row 

zofs (ofs) “sumsq/ (60* ( 150-ofs) ) 
next  offset 

The  minima  among  all  the  zofs  are  the  layers.  They  are  in  order  from 
highest  altitude  to  lowest  (since  ofs  increases  with  decreasing 
altitude).  Those  ofs  which  are  very  large  (a  low  layer)  are  cast  out, 
since  the  number  of  point*  is  insufficient  to  be  reliable,  and  the  best 
five  of  those  remaining  are  taken  to  be  the  layers.  There  may  be  only 
one  or  evon  zero  layers,  if  the  sky  is  clear. 

The  reason  the  minima  are  the  layers  can  be  seen  from  a  considerat ion 
of  either  array.  Each  is,  if  the  images  did  not  overlap,  a  purely 
random  collection  of  points  relative  to  the  other,  so  that  at  any 
arbitrary  shift  (the  value  of  ofs),  the  difference  in  magnitude  is  a 
normal  random  variable.  Hence  the  zofa  for  that  shift  will  be  close  to 
the  mean  of  the  array. 

But  at  some  offset,  a  number  of  points  will  line  up|  their  difference* 
are  no  longer  random.  Hence  zofs  for  that  offset  is  biased  downward. 
Since  clouds  are  in  general  bands  of  points  or  roughly  the  same 
magnitude,  a  plot  of  zofs  with  ofs  is  largely  a  smooth  curve  but  with 
sharp  bottoms.  The  offsets  at  which  these  bottoms  occur  are  the 
locations  of  layers  and  replace  the  term  (U1-U0)  in  (1)  above.  That  is 
(U1-U0) *minofs  where  minofs  is  the  value  of  k  at  a  local  minimum  of 
zofs (k) ,  k«0. .150. 

The  lowest  layer  that  can  be  corresponded  reliably  can  be  computed  from 
(1)  above  as  720#Sx/280  or  2.6  times  the  distance  between  the  cameras. 
This  altitude,  2600  feet,  is  too  large  for  judging  the  height  of  low 
ceilings  as,  for  example,  in  airport  operations,  so  that  the  cameras 
need  to  be  brought  closer  together  and  the  focal  length  of  the  lens 
made  shorter. 

Both  these  approaches  have  negative  aspects.  Shortening  the  base  line 
reduces  the  accuracy  in  measuring  the  height  of  very  high  clouds,  while 
decreasing  focal  length  (the  practical  optical  limit  due  to  reflection 
is  now  8  mm)  gains  only  a  factor  of  2.  ft  combination  of  an  8  mm  lens, 
and  a  base  line  length  of  155  feet  gives  a  minimum  altitude  of  200 
feet,  quite  usable.  Unfortunate) y,  at  23,000  feet,  the  one-pixel 
resolution  is  about  50*  of  the  altitude. 

ft  solution  for  low  altitudes  which  has  not  been  worked  out  is  a  time 


(2) 

(3) 

(4) 
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series  approach.  In  this  approach,  a  picturs  is  taksn  from  on*  camera 
and  th*  radial  velocity  of  th*  cloud  lay*r  determined.  In  principl*, 
th*  pictur*  will  not  changa  greatly  as  it  moves  across  th*  sky,  and 
eventually  some  part  of  it  will  appear  in  view  of  the  other  camera 
(unless  it  is  moving  directly  across  the  base  line.)  Then,  being  seen 
again,  some  estimate  of  its  altitude  can  be  mad*  by  a  sort  of  time 
laps*  correspondence  process. 

A  fourth  option  and  probably  the  best  compromise  is  to  maintain'  a 
baseline  sufficient  to  give  required  accuracy  at  high  altitudes,  then 
tilt  the  cameras  toward  each  other  until  reliable  overlap  occurs  for  a 
layer  at  2001 .  Distance  apart  (Sx)  for  overlap  at  such  an  altitude 
depends  on  focal  length,  FL,  number  of  pixels  in  th*  x  direction  of  the 
chip,  Npx  and  angle  of  tilt,  Ai 

Sx(Zmin-200» )»2#Zmin/tan(90-2*atn(Npx/2/FL> 

where  tilt  angle  A*Npx/2/FL  equals  half  the  apex  angle  so  that 
the  zenith  is  visible  at  on*  end  of  th*  chip. 

Accuracy  of  an  observation  of  altitude  is  the  fraction  *  where 

e-(Z(du)-Z(du+l) )/Z(du) 

and  du  is  th*  offset.  Here  always,  du«(l~*)/e.  Then  the  altitude 
above  which  the  accuracy  of  the  reading  is  less  than  *  is 

Z(*>»(FL-Sx>**/(l-*> 

Since  th*  field  subtended  by  Npx  pixels  decreases  as  FL  increases, 
altitude  for  a  given  error  and  given  Zmin  does  not  decrease  a*  fast  as 
FL  ncreases.  Doublinq  focal  length,  from  B  to  16  m  say,  halves  Sx 
from  533*  to  213',  but  decreases  the  5*  Z  from  10103*’  to  8083*,  only 
20*.  Further  as  focal  length  increases,  spherical,  aberration 
decreases,  and  more  of  the  chip  can  be  used  without  compensation. 

Hence  for  low  altitude  work,  a  reasonable  siting  policy  is  to  select  a 
longer  focal  length,  a  shorter  base  line  and  a  tilt  equal  to  half  apex 
angle  and  accept  imprecision  at  higher  altitudes.  If  high  altitude 
accuracy  is  required,  a  third  camera,  C,  is  added  to  th*  end  of  the 
line  A,  3,  C  located  so  that  its  field  of  overlap  with  A  begins  where  th* 
accuracy  of  the  pair  AB  drops  below  the  acceptable  level.  That  is.  its 
location  is 

Sx  (AC)  «C2*  <*/  ( 1-e)  )  *FL*2*Zmin3  /  Ctan<^0-2*atn(Npx/2/FL>  **23 

Th*  necessary  transformat ion  of  coordinates  of  a  point  U,  V  in  th* 
tilted  chip  plane  to  a  point  Q,  R  in  an  imaginary  horizontal  plan*  is 
computed  as  follows! 

h~  FL/cos(L)  where  A  is  th*  tilt  angle,  and  FL  the  /oca!  length, 

B"  atn<U/FL)  where  U  is  th*  x  coordinate  of  th*  point  in  th*  chip. 

Q»  h*tan(A-B)  where  Q  is  the  new  x  coordinate  of  the  point  and 
R»  V*h/ <FL*cos  (A-B)  /cos  (B) ,  where  R  is  the  new  y  coordinate. 
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The  overlaying  process  is  applicable  in  the  calculation  of  dx  and  dy, 
the  velocities  of  the  layers.  First  however,  the  points  that  "lined 
up"  when  a  given  layer  was  found  must  be  identified.  These  are,  of 
course,  those  points  with  small  del,  since  these  small  dels  were  what 
reduced  zof  s  in  the  first  place.  Hence,  one  zimg  array  is  placed  at 
that  layer's  ofs  over  the  array  from  the  other  camera  for  the  same  time 
and  all  points  for  which  del  is  "close"  to  zero  are  copied  into  a  new 
array,  vimg<0, i,j),  while  the  others  are  set  to  the  mean  for  the  array. 
"Close"  is  taken  arbitrarily  to  be  the  square  root  of  variance  of  the 
values  of  the  array  around  its  mean,  about  one  standard  deviation. 

If  a  copy  occurs  and  the  value  is  not  "black"  (meaning  blue  sky  can  be 
seen),  a  counter  "cover”  is  incremented,  and  eventually  divided  by  the 
number  of  points  compared,  giving  the  fraction  of  the  sky  that  contains 
a  cloud  at  this  level  provided  it  is  not  obscured  by  a  cloud  at  a  lower 
level.  Later,  in  the  Report  writing  routine,  the  fractions  are  changed 
to  the  standard  meaning  of  cover. 

The  layer-f inding  procedure  is  performed  on  the  arrays  zimg  from  both 
cameras  for  the  next  time  period,  producing  another  array  vimg(l,i,j), 
where  again,  the  first  index  is  the  time  period.  Now  vimg(0)  is 
shifted  over  vimg(l),  in  a  rectangular  outward  moving  spiral  pattern 
until  a  minimum  is  found  that  is  not  replaced  for  two  full  circles. 
The  assumption  is  that  there  is  only  one  minimum  dx,dy  in  the  layer, 
but  if  there  are  more  than  one,  the  true  one  is  likely  to  be  the  one 
nearest  0,0.  The  algorithm  is  described  in.  more  detail  in  Annex 
3. 2.  8.  2. 

The  terminal  values  of  du  and  dv  are  then  passed  along  to  Report () 
where  they  are  adjusted  for  pixel  width  and  altitude  of  the  layer  to  dx 
and  dy,  and  velocity  and  direction  are  computed  approximately  as 

velocity (layer)  ■  sqrt (dx*dx+k*dy*dy> | 

direct  ion (layer)  ■  atn(dy/dx)+angle  of  the  base  line. 

(Accommodation  is  made  for  the  pair  of  equal  values  of  atn  in  the 
circle  and  for  the  fact  that  dx  or  dy  may  be  zero. ) 

An  error  exists  in  these  algorithms  for  the  following  reason.  If  the 
lenses  were  pinholes  and  the  pixels  square,  the  points  that  line  up 
would  have  displacements  from  the  center  of  the  camera  exactly 
proportional  to  their  displacements  in  x  and  y  from  the  "master" 
camera.  The  non-squareness  of  the  pixels  (1.17  to  1)  is  easily  fixed, 
but  the  lens  is  not  a  pinhole.  Hence  motion  of  points  in  the  periphery 
of  the  image  for  unit  velocity,  which  is  assumed  constant  for  all 
points  in  the  sky  at  the  same  altitude,  is  less  than  the  motion  of 
points  in  the  center  of  the  image.  Hence  there  is  some  "blurring"  of 
the  overlay,  since  the  entire  image  is  compared.  However,  the  error 
appears  to  be  small — velocities  are  perhaps  a  percent  low,  and 
altitudes  about  a  percent  high. 

There  does  exist  a  "worst  case"  condition  under  which  the  dxdy 
algorithm  will  fail.  If  the  layer  consists  of  one  small  cloud  moving 
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at  more  than  its  radius  in  two  seconds,  there  will  be  no 
correspondence.  However,  in  this  case,  either  the  cloud  is  so  small  as 
to  be  insignif icant,  or  else  it  is  very  low  and  moving  so  fast  that  the 
instruments  themselves  will  likely  have  blown  away.  This  does  remove 
the  problem  of  landing  aircraft,  birds,  bugs  and  blowing  leaves, 
however.  If  the  problem  is  significant,  the  inter-sample  interval  can 
be  reduced. 

Percent  cover  is  simply  100  time?  the  ratio  of  number  of  cells  in 
velocity  image  which  are  not  "black"  to  the  total  number  of  cells, 
9000,  in  the  array.  For  each  layer,  the  fraction  of  cover  is  1  minus 
the  ratio  of  the  number  of  cells  not  included  in  this  layer  but 
included  in  lower  layers.  An  error  that  exists  here  is  that  the  image 
is  not  of  the  w^ole  sky,  but  rather  of  a  rectangular  cone  about  28 
degrees  in  apex  angle.  Over  time,  however,  if  there  is  any, wind  aloft, 
clouds  in  an  unobserved  part  of  the  sky  will  move  over  the  instrument 
and  become  visible,  so  that  a  time  integral  becomes  a  suitable  measure. 

The  computation  of  standard  cover  proceeds  as  followsi 

cover (this  layer)  ■  cover (next  lower  layer) 

cover (this  layer)* (1-cover (next  lower  layer)) 

where  cover  is  the  fraction  derived  in  the  function  "makevimg" 
described  above.  There  is  some  question  that  this  is  in  fact  a  good 
estimate  of  cover,  particularly  for  high  scattered  layers  above  broken 
layers.  More  experience  with  skies  with  .  these  particular 
characteristics  will  be  needed. 


4  History  of  the  Phase  I  Study 

The  work  of  the  Phase  I  study  was  performed  in  four  stages,  essentially 
as  described  in  the  proposal. 

The  specifications  of  a  number  of  electronic  cameras  were  evaluated  and 
the  choice  was  made  to  employ  the  General  Electric  TN2305  CIO  units 
which  seemed  to  have  the  best  character ist ics  and  the  highest 
evaluation  by  users  for  sensitivity,  uniformity,  range  and  absence  of 
artifacts  in  the  images  produced.  A  system  integration  firm,  Poynting 
Products  Inc,  of  Chicago  packages  the  grabber  boards  made  by  Epix,  Inc 
with  their  software  and  GE  cameras  for  the  surveillance  market.  After 
some  negotiation  and  preparation  of  specifications  for  the  WHSKY 
application,  two  units  were  ordered  subject  to  the  ability  of  the 
cameras  to  transmit  the  data  600  feet  in  coaxial  cable.  Experiments  by 
Poynting  beinp  favorable,  the  equipment  arrived  at  Micro  Science  in 
early  October. 

Meanwhile,  a  computer,  the  PCs  Ltd  206-12,  was  purchased  by  mail  order. 
After  delivery  it  was  assembled  and  the  various  software  system 
packages  installed.  A  "C"  interpreter  and  Lattice  C  compiler  were 
procured  and  installed.  The  principal  investigator  then  integrated  the 
Run  C  Professional,  Lattice  C,  SVPC1P  (Epix’s  demonstrat ion  program) 
with  the  computer,  arranged  for  and  supervised  the  survey  of  the  base 
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line,  constructed  the  stands  and  fences  around  them,  procured,  laid  and 
buried  4000  feet  of  coaxial  and  110  v.  power  cable,  and  designed  the 
Mounts  for  the  cameras. 

It  was  originally  anticipated  that  the  mounts  would  be  constructed  of 
jig  stock  by  a  machine  shop,  but  the  cost  was  unacceptably  high,  so  the 
PI  fabricated  tnerr.  in  Micro  Science’s  shop  using  high-quality  plywood. 
Os  it  turned  out,  the  precision  and  rigidity  of  the  mounts  went  more 
than  sufficient  for  the  current  effort.  Intermittent ly  during  this 
period,  the  algorithms  of  the  simulations  were  recoded. 

When  the  Poynting  hardware  and  system  software  arrived,  the 
boards  were  installed  in  the  computer  and  connections  made. 
Unfortunately  one  plug  was  inserted  upside  down  blowing  one  of  the 
boards  which  had  to  be  replaced.  Later,  a  different  kind  of  failure  on 
another  board  required  it  to  be  replaced.  Then  at  various  timns  both 
cameras  had  to  be  replaced. 

P  very  small  portion  of  the  Epix  software,  called  SVPC1 P,  is 
actually  relevant  to  the  WHSKY  project.  Pll  of  the  functions  used 
except  invocation  of  the  grabber  boards  themselves  have  been  replaced 
by  Micro  Science  programs}  the  major  effort  has  been  to  figure  out  how 
to  get  at  the  boards  without  having  access  to  their  source  code,  which 
is  for  sale,  but  at  a  prohibitively  high  price.  If  Poynting  ar»d  Epix 
are  to  be  used  again,  particularly  for  a  commercial  product,  a 
substantial  negotiation  will  be  required  to  decide  how  to  split  down 
the  interests. 

During  October  and  early  November,  attempts  were  made,  first 
with  one  board  then  the  other  and  one  camera  then  the  other  to  solve 
two  apparently  nearly  intractable  problems.  One  was  to  adjunt  the 
image  in  Such  a  way  that  the  sky  as  seen  through  a  standard  lens  (which 
simulates  the  eyeball)  was  transformed  to  appear  as  if  it  were  seen 
through  a  pinhole.  The  reason  for  doing  this  is  discussed  in  detail  in 
Pnnex  3.4.2  Pdjust lens ( ) .  The  second  problem  was  to  find  a  way  to 
handle  substantial  variation  in  the  signals  from  the  same  sky  its  seen 
by  two  cameras  1000  feet  apart.  The  solution  is  also  described  in 
Pnnex  3.4.1.  FiiterO.  Ps  it  has  turned  out,  neither  solution  was 
required,  since  the  problems  are  insignif icant,  easily  solved  in  much 
simpler  fashion. 

By  mid-November,  the  instrument  itself  consisting  of  the  two 
cameras,  the  mounts,  the  stands,  the  cables,  the  boards  and  the 
computer  was  completely  assembled,  installed  on  the  test  bed  and  made 
operational.  Serious  development  of  the  algorithms  began.  Pjthough 
the  development  went  through  some  20  editions  of  the  code,  th«>  final 
product,  which  is  described  in  Pnnex  III  is  essentially  whut  was 
planned  in  November.  The  demonstrat ion  program  presented  iH  early 
January  made  liberal  use  of  "filter".  This  has  since  been  scrapped  and 
replaced  by  modifications  of  the  two  major  overlay  programs  which  have 
the  effect  of  producing  filtered  images.  It  is  possible  that  in  some 
cases  of  layers  with  relatively  few  points  (the  layer  is  Heavily 
obscured  by  lower  layers)  that  some  kind  of  filtering  routine  will  be 
needed  to  reduce  the  number  of  minima  (which  give  dx  or  dy)  to  one.  To 
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date,  no  cloud  pattern  has  been  recorded  which  is  not  tractable  to 
existing  algorithms. 


By  the  end  of  January  the  program  was  complete  except  for 
cleanup  and  rr<*organizat ion.  The  current  operational  version  is  listed 
and  discussed  in  Annex  3.  Due  to  the  delays  and  problems  described 
above,  plus  an  unanticipated  paucity  of  the  kind  of  skies  which 
contained  the  clouds  of  interest,  extremely  cold  weather  in  December 
and  two  feet  of  snow  in  January,  studies  of  summer  cumulus, 
thunderstorms,  lenticulars  and  other  such  exotic  forms  did  not  take 
place. 


On  the  other  hand,  a  great  deal  of  work  was  done  with  stratus 
and,  slant-range  from  indoors,  on  snowing  clouds.  The  cameras  do  not 
work  well  below  freezing,  and  since  they  are  not  sheltered  cannot  be 
left  but.  The  kinds  of  problems  which  remain  to  be  dealt  with,  as 
described  in  the  proposal  for  stage  3  of  this  Phase  I,  are  addressed  in 
the  new  proposal  for  Phase  II. 


In  summary,  the  project  did  not  present  any  particularly 
unexpected  or  unsolvable  problems,  and  has  been  completed  within  the 
original  time  period,  within  funds,  even  though  a  no-cost  extension, 
requested  when  such  an  outcome  did  not  seem  too  likely,  was  granted. 


S  Findings  of  the  Phase  I  Study 

While  the  primary  goal  of  Phase  I  was  to  get  an  instrument 
assembled  and  into  operation  with  existing  algorithms,  which  was 
achieved,  a  number  of  discoveries  were  made  which  are  of  interest. 
These  were  primarily  of  the  form,  in  what  way  does  the  real  sky  differ 
from  the  simulated  skies  that  have  been  used  heretofore?  The  general 
answer  is,  in  quite  a  few  ways,  but  in  the  direction  to  make  analysis 
easier  rather  than  harder  than  expected. 

Admittedly,  part  of  this  simplification  was  the  result  of 
certain  unexpected  characterist ics  of  the  overlay  algorithms  whose 
behavior  has  only  recently  been  understood.  The  complexity  of  the 
problem  turns  out  to  be  smaller  than  anticipated  as  well  because  of  the 
number  of  points  in  the  sky  directly  overhead,  that  is,  in  the  cone  of 
view  <28  deg.  apex  angle  around  the  zenith). 

Since  we  are  not  normally  in  the  habit  of  looking  straight  up 
at  the  sky,  we  imagine  cloud  structures  to  be  complex,  as  seen  in  the 
first  20  or  so  degrees  above  the  horizon  where  we  see  the  sides  of 
clouds  and  a  far  greater  expanse  of  sky  and  number  of  clouds  than  is 
actually  directly  overhead.  Hence  a  given  observation  is  a  very  small 
sample  of  the  visible  sky,  and  therefore  quite  homogeneous. 

Subsequent  studies  will  want  to  combine  time  series  of  the 
local  area  of  visibility,  developing  changing  patterns  from  a  series  of 
individual  reports  rather  than  trying  to  grasp  the  entire  sky  from  one 
picture. 
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Another  discovery  of  interest  is  that  because  the  CID  chips  in 
the  cameras  have  a  much  wider  spectrum  of  frequency  sensitivity,  cloud 
layers,  particularly  stratus,  which  appear  completely  blank  to  the 
naked  eye,  do  in  fact  have  quite  elaborate  structure,  allowing  for 
triangulat ion  under  unexpected  conditions.  These  may  be  holes  through 
whic*~>  the  cameras  see  overlying  layers.  Even  a  completely  clear  sky 
will,  at  sunset  or  dawn,  display  a  granular  image,  this  being  the 
reflection  of  sunlight  from  the  ridges  and  surfaces  of  the  top  of  the 
smog  layer  which  is  normally  invisible. 

Furthermore,  clouds  can  be  seen  even  when  light  snow  or  rain  is 
falling,  and  experiments  with  the  time-exposure  feature  of  the  cameras 
show  that  stars  themselves  are  visible  on  a  dark  clear  night.  Pixel 
electronic  noise  is  severe,  however,  and  further  work  needs  to  be  done 
to  implement  methods  to  cancel  out  the  noise. 

It  is  certainly  clear  that  cameras  must  be  equipped  with  lenses 
with  auto-iris.  The  cameras  and  the  Poynting  hardware  make  provision 
for  the  control  of  "motorized"  irises.  It  is  simply  a  matter  of  adding 
the  feature. 

The  cameras  and  their  stands  must  be  sheltered  in  some  way.  We 
propose  a  shelter  resembling  a  miniaturized  version  of  the  conventional 
Post  Office  street  corner  mailbox,  which  has  a  pull-down  letter  chute 
and  a  round  top.  Suppose  the  door  were  hinged  to  the  top  of  the  slot 
and  the  camera  mounted  inside  looking  nov  straight  up,  but  at  a  small 
angle  through  the  doorway.  When  it  was  time  to  take  a  picture,  a  fan 
inside  would  start,  blow  the  door  open  and  the  airstream  would  blow 
away  any  snow  or  rain  or  other  foreign  matter  from  the  opening.  After 
the  picture  was  taken,  the  fan  would  turn  off  letting  the  door  close. 
Some  arithmetic  (see  2.6)  would  be  required  to  perform  the  necessary 
transformat  ion  of  the  image  into  what  would  be  seen  if  the  axis  were 
perpendicular,  but  other  than  that  there  seem  to  be  no  problems  of 
significance. 

The  actual  application  and  system  software  package  for  an 
operating  system  can  be  much  smaller  than  that  used  in  the  development 
system.  No  hard  disk  would  be  required,  for  example,  but  rather  the 
code  could  be  written  to  ROM  (read-only  memory)  and  a  much  simpler 
computer  used.  Likewise,  there  is  no  need  for  2  megabytes  of  grabber 
ram,  nor  a  monitor.  In  the  most  likely  application,  the  signals  would 
be  transferred  by  telephone  line  to  the  point  of  use,  or  combined  with 
other  messages  from  an  automatic  weather  observation  station.  A  large 
block  of  ram  could  be  provided  in  the  hardware  to  accumulate  data  from 
successive  pictures,  more  or  less  like  a  flight  recorder,  so  that  the 
newest  data  simply  overwrites  the  oldest.  This  would  allow  rather 
sophisticated  time  series  analyses  to  be  made  without  requiring  any 
mechanical  storage  device. 

Finally,  we  sense  from  the  kinds  of  problems  encountered  that 
the  algorithms  are  quite  robust,  that  is,  not  easily  thrown  off  track 
by  unexpected  conditions.  A  great  deal  of  data  is  thrown  away  or 
aggregated,  but  the  sky  is,  at  least  to  the  degree  that  it  has  been 
studied,  not  a  difficult  phenomenon.  While  it  is  a  statistical 
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phenomenon,  in  aggregate  the  characterist ics  of  interest  are  easily 
found.  In  spite  of  their  apparent  chaotic  and  turbulent  nature,  clouds 
are  in  fact  computationally  tractable. 


&  Conclusions 

P  prototype  Whole  Sky  Sensor  has  been  built  and  tested  and 
while  there  has  been  a  very  limited  amount  of  real  time  utilization, 
has  proven  adequate  for  the  task  of  determining  height  and  velocity  of 
layers  of  clouds.  Considerable  additional  research,  particularly  into 
cloud  types  and  formations  of  interest  to  the  micro-meteorologist  needs 
to  be  done,  and  the  device  needs  a  suitable  shelter. 

The  next  stage  of  development  should  be  directed  toward 
answering  these  requirements,  with  extensive  field  testing  under  user 
management,  further  algorithmic  development,  and  interfacing  with 
existing  met  systems.  Ps  a  part  of  this  phase,  it  should  be  possible 
to  design  a  unit  for  general  commercial  production  and  instal lat ion. 
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Annexes 


Alt  Equipment  specifications 
Al. 1  Cameras 

The  General  Electric  4TN2505A3  Solid  State  Surveillance  cameras  contain 
a  Charge  Injection  Device  (CID)  chip  with  a  248(V)  by  388(H)  pixel 
array.  The  pixels  themselves  are  27.3  uM(V)  by  23.3  uM(H).  The  width 
of  the  usable  portion  of  the  chip  (360  pixels)  is  8.388  mm.  Thie 
divided  by  the  focal  length  of  a  16  mm  lens  gives  the  subtended  angle 
to  be  arctan (8. 388/16) *27. 66  degrees.  At  10,000  feet  above  the  camera, 
the  width  of  the  field  is  thus  5242  feet  long  (along  the  x-axis)  and 
240/360*1. 1 17*. 32425  *  3904  feet  wide.  A  pixel  subtends,  at  this 
height,  in  x,  5242/360  »  i5  feet,  and  in  y,  16  feet. 

Relevant  pages  from  manufacturer’ s  data  sheets  are  attached  at  the  end 
of  this  Annex. 


A1.2>  Frame  grabbers 

There  are  two  grabbers,  identical  except  that  one  is  programmed  to 
respond  to  signals  directed  to  the  "master"  and  the  other  as  "slave". 
This  discussion  describes  either  one  of  the  units,  except  for  that 
difference. 

The  grabber,  manufactured  by  Epix,  Inc.  is  a  printed  circuit  board 
which  occupies  one  of  the  two  short  slots  in  the  computer,  sharing  the 
data  bus.  It  is  powered  by  the  computer  power  supply  and  SNC 
connectors  provide  access  for  signals  entering  and  leaving  the  grabber. 
The  board  includes  a  flash  8-bit  ADC,  a  ROM  chip  holding  certain 
proprietory  software,  a  RAM  chip  by  which  instructions  for  a  specific 
application  car.  be  set  into  the  grabber  from  the  computer,  1  megabyte 
of  data  RAM,  and  the  usual  ancillary  components. 

When  the  board  is  invoked,  parameters  are  stored  in  the  RAM  and  if  set 
to  do  so,  the  analog  TV  signal,  which  streams  from  the  camera  into  the 
AEu  continually  as  long  as  the  camera  is  powered,  is  routed  as  a  string 
of  numbers  from  the  ADC  to  a  location  in  the  grabber  data  ram.  Every 
frame  (equivalent  to  a  TV  picture  frame)  the  previously  written  data  is 
overwritten  by  new  data. 

However,  i  software  function  can  be  used  to  change  the  starting  address 
of  the  frame  so  that  N  pictures  can  be  written  to  ram  where  N*  1048576 
modulo  frame  size,  in  bytes,  so  that  a  kind  of  "movie"  can  be 
accumulated,  or  several  frames  in  sequence  stored  for  later  analysis. 

The  picture  is  simultaneously  passed  through  a  digital  to  analog 
converter  (back  to  its  original  form)  and  output  «s  a  signal  available 
to  drive  a  standard  TV  set,  the  monitor,  on  which  the  user  can  observe 
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the  picture  being  taken. 

Once  the  desired  picture  is  in  the  grabber  ram,  the  input  to  the  area 
is  disabled  and  the  resulting  block  of  memory  can  be  transferred 
elsewhere  in  the  computer  or  stored  to  disk. 

A  very  large  number  of  abilities  have  been  programmed  for  this  device 
by  Epix,  but  most  are  unused  in  this  appl  icat  ion.  Relevant  data  sheets 
from  the  manufacturer’ s  document  at  ion  are  attached  at  the  end  of  this 
Annex. 


A1.3s  Computer 

The  PC’s  Ltd  286-12  computer  acquired  for  this  project  is  the  fastest 
of  their  family  of  IBM  Personal  Computer  compatibles,  and  as  of  last 
August,  the  fastest  such  machine  available  from  anyone.  By  Epix's 
algorithm  (which  determines  the  clock  and  processing  rates  of  the  CPU 
in  order  to  manage  its  own  firmware)  the  286-12  is  very  nearly  4  times 
as  fast  as  a  standard  8086-based  IBM  PC,  twice  as  fast  as  an  IBM  AT  and 
is  a  third  the  price  of  the  latter.  Nonetheless,  it  h«is  proven 
entirely  reliable,  and  the  12  megahertz  clock  rate  has  made  child’s 
play  of  the  necessary  number  crunching  and  array  manipula.t ion. 

Relevant  pages  from  the  manufacturer*  s  specifications  are  attached  at 
the  end  of  this  annex. 


A1.4i  Installation 

The  cameras  are  mounted  in  semi-gimballed  frames  which  are 
placed  on  the  stands  at  either  end  of  the  baseline.  Since  these  may 
not  be  visible  from  one  another,  adjustments  need  to  be  made  to  assure 
that  the  cameras  are  in  line  and  looking  straight  up.  This  latter  is 
done  by  "swinging"  the  cameras  on  their  gimbals  so  as  to  achieve 
registration  at  various  angles  along  the  line,  thus  assuring  there  is 
no  tilt  of  the  z-axis  around  the  x-axis,  and  that  the  x-axis  passes 
through  the  centers  of  both  cameras. 

Then  the  z-axes  are  set  perpendicular  by  means  of  a  float 
level.  Error  in  this  setting  appears  as  error  in  altitude  of  layers. 
It  also  appears  as  a  "smearing"  of  the  bottoms  of  the  overlay  display, 
since  the  effect  of  an  error  is  as  if  the  sky  were  tilted  along  the 
x-axis  in  either  or  both  of  the  cameras.  A  very  precise  alignment  can 
be  made  from  stars,  since  the  rays  from  a  star  to  each  of  the  cameras 
are  essentially  parallel.  Two  images  should  exactly  match  (within  the 
resolution  of  a  pixel)  if  the  cameras  are  set  correctly. 

A  very  precise  measurement  of  perpendicularity  of  the  z-axis 
can  be  made  by  the  method  of  bort-si ght ing.  In  this  method,  a  straight 
tube  one  inch  or  less  in  inside  diameter  10  feet  long  is  suspended 
directly  over  the  camera  by  the  top  end.  rhen  if  the  camera  is  set 
correctly,  the  image  of  the  opening  at  the  top  of  the  tube  is  visible 
on  the  monitor  through  the  bottom  of  the  cube.  Assuming  that  the 
images  can  be  located  within  one  tenth  inch,  the  angular  error  out  of 
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perpendicular ity  is  not  more  than  1  part  in  1200,  or  3  minutes  of  arc, 
which  is  less  than  the  optical  width  of  one  pixel.  380  pixels  subtend, 
with  a  18  mm  lens,  28  degrees  or  4.8  minutes  per  pixel.  Hence  the 
boresighting  method  achieves  an  accuracy  which  exceeds  the  camera’s 
resolution. 


Annex  2:  System  Software  specifications 


A2. 1  MS-DOS 

This  package,  the  Microsoft  Disk  Operating  System,  v.  3.12,  was 
purchased  with  the  computer  and  is  tailored  for  the  PCsLtd  286-12,  a 
machine  closely  related  to  the  IBM  PCAT  from  whose  specif ications  it 
has  grown.  At  least  part  of  the  virtue  of  RCP  (see  below)  is  the 
"shell"  facility  which  allows  DOS  functions  to  be  called  from 
application  programs  running  under  RCP.  Many  of  these  functions,  in 
particular  copy  and  concatenate,  used  to  combine  and  store  collections 
of  files  containing  segments  of  pictures,  get  considerable  exercize. 


A2. 2 i  Lattice-C 

The  C  compiler  language,  written  by  Kernighan  and  Ritchie  of  Bell  Labs 
in  the  late  ’70’s  primarily  for  the  development  of  the  Unix  operating 
system  has  become  the  language  of  choice  for  system  developors,  quickly 
superceding  Pascal,  PL/1,  Fortran,  APL,  BASIC  and  assembly  language  as 
the  tool  to  use  in  the  development  of  compact  highly  efficient  easily 
written  and  error-resistnat  code. 

Because  of  its  popularity,  several  C- interpreters  have  emerged  in  the 
last  few  years  to  provide  C  programmers  with  a  tool  by  which  programs 
can  be  written  and  executed  in  the  interpretive  mode,  and  then  compiled 
into  machine  language  after  all  the  "bugs”  have  been  removed. 

This  combination  of  a  powerful  language  in  both  the  compiled  and 
interpretive  forms  with  easy  links  between  the  two  and  with  a  large 
supply  of  off-the-shelf  program  packages  and  programming  aids,  viz., 
"C-food  Smorgasbord"  has  made  software  development  rapid  and  accurate, 
perrr  tting  the  analyst,  who  figures  out  what  the  computer  is  to  do,  to 
wr;.te  his  own  programs  and  debug  them,  module  by  module  as  he  generates 
the  algorithms  and  procedures  to  be  implemented,  then  compiling  the 
specific  modules  into  machine  language,  and  calling  them  from  the 
interpreter  as  required. 

This  procedure  is  very  fast  and  very  efficient,  and  has  been  employed 
in  this  development.  If  there  is  any  problem  it  is  that  programs  are 
so  easy  to  write  and  test  that  they  pile  up,  and  keeping  track  of  the 
procedures  that  are  available  becomes  a  bit  of  a  bookkeeping  problem, 
since  large  harddisk  memories  tempt  the  programmer  to  save  functions 
which  ought  to  be  discarded,  simply  on  the  grounds  that  they  might 
sometime  be  useful.  But  this  is  a  small  price  to  pay  for  the  fact  that 
rather  large  programs  can  be  turned  out  in  short  order  by  a  single 
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analyst  working  alone,  once  he  has  mastered  the  language. 

There  are  a  number  of  texts  on  "C"  which  make  the  transition  to 
the  language  not  necessarily  easy  but  achievable  by  anyone  with  a 
reasonable  amount  of  programming  experience,  even  if  it  has  been 
limited  to  BASIC,  an  interpretive  language.  I  strongly  recommend  that 
any  scientist  or  engineer  who  has  a  need  for  computation  try  this 
package 

The  Lattice  C  compiler  is  neither  the  easiest  nor  the  fastest  C 
compiler  available,  but  it  is  probably  the  most  popular,  in  terms  of 
the  number  of  lines  of  code  that  have  been  written  in  it.  .Lattice  was 
used  to  write  SVPCIP  and  Run-C  Professional  with  which  it  interfaces 
nicely,  and  the  large  memory  and  -k2  options  are  available  for  the 
Intel  80286/80287  processor  chip  pair  which  exist  in  the  PC's  Ltd 
computer. 


A2.3:  Run/C  Professional 

This  C  interpreter,  written  by  Walton  and  Brooks  of  Lifeboat  Associates 
is  very  nearly  flawless  as  a  program.  At  every  point  at  which  I 
imagined  that  the  interpreter  had  made  an  error,  it  turned  out  to  be  my 
programming  error.  If  there  is  a  criticism  it  is  that  certain  kinds  of 
errors,  which  should  be  reported  as  syntax  errors,  produce  major 
crashes  or  completely  unrelated  diagnostics. 

Program  size,  particularly  data  arrays  of  the  sort  used  in  WHSKY,  can 
be  cramping,  in  part  because  a  third  of  the  memory  is  blocked  off  from 
use  by  the  grabber  boards,  and  by  the  way  RCP  organizes  memory.  On  the 
other  hand,  with  practice  it  is  almost  always  possible  to  write 
programs  which  run  the  first  time  and  their  compilation  takes  a  couple 
of  minutes  at  most.  This  allows  the  prgrammer  to  concentrate  on  the 
problem  to  be  solved  rather  than  the  means  of  solving  it,  *  highly 
desirable  characterist ic  of  any  language. 

A2. 4 i  SVPCIP 

This  software  package,  written  in  C  by  Epix  Associates  is  designed 
primarily  to  be  used  in  testing  and  evaluating  the  Silicon  Video  frame 
grabber,  and  to  serve  as  a  basis  for  customizing  applications  in  the 
surveillance  area.  It  is  not  a  stand-alone  program,  requiring  that  the 
user  enter  menu  selections  from  the  keyboard.  It  contains  a  large 
number  of  functions  of  various  sorts,  such  as  different  kinds  of 
digital  filters  to  be  applied  to  the  pictures  to  enhance  viewing,  for 
changing  frame  size,  for  rearranging  parts  of  the  pictures,  zooming  and 
the  like. 

The  package  as  a  whole,  albeit  a  very  respectable  bit  of  work,  is 
overkill  for  this  project,  but  still  prevents  stand-alone  operation 
which  ought  to  begin  automatically  with  power  on,  and  continue  without 
intervention  until  power  down.  Nonetheless,  SVPCIP  has  served  this 
Phase  I  effort  well  and  without  question  the  minor  inconveniences  can 
easily  be  eliminated,  at  a  price. 


m 
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Annex  3*  The  WHSKY  Computer  Program 


A3. 1  Overall  flows 

In  the  normal  mode  of  use,  the  instrument  runs  continuously,  making 
observations  every  ten  minutes  as  long  as  power  is  on.  But  to  get  the 
system  going,  certain  steps  are  required  of  the  operator.  First  both 
cameras  are  turned  on,  and  the  computer  powered  up.  DOS  boots  control 
directly  to  the  RCP  interpreter  by  an  Autoexec.bat  invocation,  RCP\RCP 
(with  sizing  parameters).  Then  operator  loads  and  runs  RCPNWHSKY. 200, 
the  current  version  of  the  system.  Whsky  wakes  up  the  grabber  software 
package,  SVPCIP,  as  a  "shell"  activity  from  within  RCP.  Control  is 
returned  to  RCP  by  the  user  entering  "qu"  and  2  carriage  returns. 

Some  further  initializing  occurs  in  this  version  before  the  main  loop 
begins.  RCP  loads  the  library  of  compiled  Whsky  modules,  called 
FEBLIB.C  and  listed  elsewhere,  MATHLIB,  a  package  of  mathematical 
functions  and  GRAPHLIB,  a  package  of  graphics  functions.  Both  MATHLIB 
and  GRAPHLIB  were  provided  with  RCP.  Then  the  set  of  parameters 
required  by  the  grabber  boards  is  copied  from  a  disk  file  to  the 
grabbers.  These  parameters  were  originally  input  manually  (from  the 
keyboard)  and  copied  to  hard  disk  by  a  small  service  routine  written 
for  the  purpose. 

Then  the  program  sends  an  eight-bit  string  of  flags,  hex  number  0xac, 
via  output  port  0x301  to,  the  grabber  boards,  telling  them  to  digitize 
the  incoming  camera  images.  The  images  also  appear  on  the  monitor,  so 
that  operator  can  adjust  the  irises  of  the  cameras.  Focus  is  set  to 
infinity.  The  program  waits  for  a  keyboard  input  signalling  thatv.  the 
cameras  are  set  to  the  correct  exposure.  Once  that  key  stroke  has  been 
entered,  the  system  runs  unattended,  looping  tnrough  the  balance  of  the 
program  unattended  as  long  as  power  is  on. 

The  first  module,  GrabimageO,  reads  a  sky  from  each  camera  at  some 
time  t0  into  the  two  grabber  rams  simultaneously,  and  does  it  again  at 
time  tl,  2  seconds  later.  Before  returning,  the  time  of  the  next 
observation,  t3,  is  calculated  and  stored.  Then  the  four  images  are 
copied  from  grabber  ram  into  computer  ram  by  module  GratozimgO,  21600 
bytes  at  a  time  and  a  subset  of  these  data  points  transferred  (and 
transformed)  to  arrays  zimg0  and  zirngl  by  module  CopyimageO. 

Then  the  routine  Zolay0<)  runs.  The  image  (now  in  zimg0C01 Cil Cjl  ) 
from  the  master  board  for  time  t0  is  overlaid  onto  the  image  (now  in 
zimg0C13 C i 3 C j 3  from  the  slave  board  for  time  t0  to  find  the  values  of 
the  difference  of  the  corresponding  pixels  for  each  offset.  The 
results  are  stored  in  an  array  zofs0Ck3,  where  k*»0. .  149.  The  process 
is  repeated  for  tl,  using  ZolaylO  and  arrays  zimglC03  and  zimglCll. 
Then  the  module  CalczO  executes,  finding  the  minima  in  the  two  zofs 
and  after  some  manipulat ion,  storing  them  in  the  two  arrays  dwCllCtl 
amd  zdwCUCtl,  where  dw  contains  the  index  (offset)  of  each  layer,  and 
zdw  the  value  of  zofs  at  that  offset.  Here,  1  is  the  layer  number,  and 
the  second  index  is  the  time  (sample  number). 
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After  the  layers  are  found,  DxdylO  executes,  taking  one  layer  at  a 
time.  First  it  makes  two  new  arrays  for  each  layer,  using  MakevimgO. 
Each  array  is  the  same  points  from  the  same  layer  from  the  master 
camera  for  two  different  times.  One  of  these  is  overlaid  in  a 
rectangular  pattern  onto  the  other  finding  the  minimum,  vofsCdul CdvD 
for  that  layer  which  represents  its  shift  in  x  and  y  in  the  time 
interval,  hence  its  velocity.  Then  du  and  dv  for  this  layer  are  stored 
in  lists  duCl]  and  dvC13  for  use  by  the  final  module  in  the  loop, 
Report ( ) . 

Report ()  spends  most  of  its  time  putting  together  an  impressive  looking 
display  on  the  screen.  It  prints  the  kind  of  report  (WHSKY),  the 
facility  location,  its  latitude  and  longitude,  height  above  sea-level, 
the  time  and  date,  and  a  column  stub  identifying  the  data  that  will 
subsequently  be  listed.  The  items  are  "cloud  base  feet  ASL",  "moving 
toward  degrees  true",  "velocity  knots"  and  "percent  cover".  Then  the 
conversions  of  the  dw,  du  and  dv  found  in  earlier  programs  are  made  to 
give  usable  values  for  each  layer.  The  formulae  are  detailed  below. 

The  reason  for  the  rather  clumsy  heading  for  wind  direction  is  that  in 
common  parlance,  winds  are  described  in  terms  of  the  direction  from 
which  they  are  blowing,  yet  clouds  and  storm  systems  and  fronts  are 
described  in  terms  or  the  direction  toward  which  they  are  moving,  e.  g., 
"a  northwest  wind  is  pushing  the  clouds  southeast",  which  can  be 
confusing. 

The  data  in  this  report  can  be  stored  as  a  set  of  numbers  rather  than 
being  printed,  or  can  be  transferred  to  a  modem  for  transmission  to  a 
manned  weather  station  for  interpretat ion  or  other  use.  The  present 
system  does  not  provide  this  nicety. 

After  the  report  is  presented,  the  program  reads  the  time  from  the 
clock  in  the  computer,  compares  it  with  time  t3,  the  next  observation, 
and  if  time  now  is  less  than  t3,  repeats  the  reading.  Eventually,  time 
now  is  equal  to  t3,  control  returns  to  the  beginning  of  the  loop  and 
GrabimageO  runs  again.  The  interval  is  usually  ten  minutes,  somewhat 
longer  than  it  takes  to  run  through  the  program. 

A3. 2  Detailed  description  of  the  program! 

It  is  recommended  that  the  listing  of  the  program,  Annex  4,  be  studied 
in  parallel  with  this  text,  Since  the  purpore  of  the  text  is  to  explain 
how  the  code  works,  and  the  how  and  why  of  certain  algorithms  which  are 
not  always  self-evident  from  the  code. 

A3. 2. 1  Preprocessing j 

When  Whsky  is  loaded  to  be  executed  under  the  control  of  RCP,  three 
options  are  implemented.  The  group  of  "#dsfine"  statements  instruct 
RCP  to  replace  each  occurrence  in  the  program  of  the  first  string  in 
the  line  with  the  second.  This  greatly  simplifies  changes  in 
parameters.  For  example,  the  maximum  number  of  layers  (nl)  that  can  be 
found  i3  currently  set  to  6,  but  it  has  been  at  various  times  in  the 
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development  1,2  or  3  as  well.  Since  nl  shows  up  many  times  in  the 
code,  it  is  useful  to  be  able  to  change  it  everywhere  with  just  one 
edit  operation. 

Npin  (number  of  pixels  in)  is  one  fourth  of  the  TV  image  of  86400 
bytes,  that  is  240  rows  by  380  columns.  Since  88400  is  a  data  object 
too  large  to  be  handled  routinely  by  any  computer  with  a  18-bit  address 
bus,  and  2**16=,65536  does  not  divide  nicely  into  86400,  and  86400/2  = 

4320O  is  too  large  to  be  expressed  as  a  16  bit  signed  integer,  and 
because  using  only  the  middle  third  of  the  image  throws  away  too  much 
usable  data  at  this  point,  the  standard  block  of  data  was  selected  to 
be  of  size  21600. 

Nsamples,  the  number  of  pictures  to  be  taken  by  a  camera  at  each 
observation,  is  actually  one  larger  than  the  number  of  pictures  used. 
The  reason  is  given  in  the  discussion  of  GrabimageO  below. 

Two, useful  ^defines  not  here  would  be  nrows  60  and  ncols  150,  since 
these  two  parameters  appear  throughout  the  program.  It  turns  out  for 
other  reasons  that  a  better  selection  of  the  dimensions  of  an  overlay 
array  might  have  been  120  by  135.  Had  these  defines  been  installed,  it 
would  be  easy  to  change  them,  and  change  them  back  if  wrong. 

Several  variables  and  arrays  are  declared  "static".  These  data  items 
are  initialized  to  zero  and  permanent  areas  are  set  aside  for  the  data 
in  memory  for  all  programs  to  use.  Arrays  or  data  items  that  are  local 
to  a  particular  subroutine  are  exactly  that.  Every  time  the  function 
is  invoked,  an  address  for  the  item  is  assigned  from  whatever  space  is 
available  at  the  time.  By  declaring  arrays  static,  their  locations  are 
permanently  fixed  and  any  function  can  get  at  them.  This  is  not  true 
of  the  compiled  functions  which  have  to  be  told  on  each  call  where  the 
data  named  in  the  function  is  stored.  , This  is  the  purpose  of  the 
string  of  "formal"  parameters  which  follows  the  name  of  a  called 
function. 

A3.  2* 2  Data  items:  (declared  static) 

time:  a  floating  point  number  containing  a  single  number  combining 
hours,  minutes,  seconds  and  hundredths  of  seconds,  generated  by  the 
function  btimerO,  which  reads  the  computer  clock  and  does  the 
arithmet ic. 

delaytime:  any  time  which  is  to  be  compared  with  the  clock. 
time2:  n^xt  observation,  usually  ten  minutes  after  this  one. 
time3:  next  sample,  usually  two  seconds  after  this  one. 
time4:  an  extended  delay  due  to  poor  seeing  (fog,  heavy  snow)  or 

transmission  difficulties.  This  would  be  set  by  "totlight'1  which  will 
be  included  when  the  cameras  are  equipped  with  autoiris,  or  by  "Report" 
if  a  modem  is  in  use. 

dt2«600:  seconds  added  to  time  now  to  get  time2. 
dt3*2:  added  to  time  now  to  get  time  3. 

dt4*1200:  twenty  minutes,  added  to  time  now  to  get  time  4. 
img0Cnpin3:  a  21600  byte  array  of  characters  (the  numbers  0. .255), 
into  which  the  grabber  panels  are  transferred  and  held  while  being 
transformed  into  the  zimg’s. 
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zimg0C23 C603 C1503,  z imgl C23 C603 Cl 503 s  Two  2  by  60  by  150  array*  of 
integers  which  hold  selected  bytes  for  each  camera  for  each  time.  The 
number  in  the  name  is  the  time,  the  first  ^.ndex  the  board,  the  second 
index  the  row  and  the  third  index  the  column.  Each  array  holds  16,000 
integers,  or  36,000  bytes.  Since  input  values  are  characters  (6  bit 
strings),  these  could  also  be  arrays  of  characters,  but  for  other 
reasons  the  bytes  coming  from  the  grabbers  are  converted  to  integers  by 
CopyzimgO  for  subsequent  use. 

vimg C23 C603 C1503 : ,  An  array  of  2  by  60  by  150  integers  which  holds 
the  values  from  zimg0C03 Ci3 Cj3  and  zimgl [03 C i 3 C j 3  that  are  identified 
as  being  erom  a  particular  layer.  The  array  is  reused  and  refilled  for 
each  layer. 

zofs0C1503,  zofslC1503:  Two  one-dimensional  arrays  holding  the 
results  of  the  z-overlaying  process,  one  for  each  time. 

vofsC3023:  a  scratch  string  used  to  hold  temporary  results  from 
the  smoothing  function. 

A3.  2.  3  MainO 

This  line  marks  the  beginning  of  the  program  proper.  Since  "C"  is  a 
structured  language,  the  Main  program,  which  must  be  included,  should 
be  short  and  limited  to  global  management  of  the  process,  which  it  does 
by  calling  functions  by  name  which  in  turn  call  other  functions,  each 
limited  to  data  to  which  it  is  given  access  by  the  caller,  or  which  is 
global  (as  described  above).  When  completely  purged  of  development 
modules,  Whsky’ s  Main  does  some  ini t ial i zat ion,  sets  up  a  loop  to  run 
forever,  and  then  within  that  loop,  calls  just  seven  functions.  These 
are  GrabimageO,  GratozimgO,  Zolay0<),  ZolaylO,  CalczO,  DxdylO  and 
Report ().  Each  of  these  will  be  discussed  in  detail  below. 

The  initialization  process  in  MainO  involves  waking  up  the  grabber 
software  by  a  "Shell"  instruction,  loading  the  libraries,  and  setting 
some  flags  to  permit  easy  exit  from  graphics  programs  (which  change  the 
screen  format)  during  development. 

Let  me  digress  briefly  here  to  present  a  lecture  on  graphics.  A 
graphics  capability,  particularly  in  systems  dealing  with  an  enormous 
amount  of  data,  is  extremely  valuable,  since  relatively  short  programs 
can  turn  thousands  of  numbers  into  simple  graphs  which  are  easy  to 
understand  and  which  display  trends  and  relationships  in  the  data  that 
are  invisible  in  the  numbers  themselves.  Another  advantage  is  that  the 
display  shows  program  errors  and  d iscont inuit ies  in  the  data,  or 
hardware  error  in  a  way  that  is  rarely  obvious  in  the  numbers 
themselves.  It  is  also  a  first  rate  marketing  or  training  tool  if 
properly  used. 

In  the  operating  program,  there  is  currently  no  requirement  for 
graphics  since  there  will  be  no  screen,  but  for  development  it  is 
almost  mandatory,  at  least  if  the  process  is  to  be  efficient. 
Regrettably  the  IBM  graphics  package,  which  is  implemented  in  the  PC’s 
Ltd  machines,  is  primitive  relative  to  the  very  classy  graphics  of  the 
Texas  Instruments  Professional  Computer,  a  very  sophisticated  cousin  of 
the  original  IBM  PC  which  did  not,  regrettably  gain  a  sufficient  share 
of  the  market  to  survive.  End  of  lecture. 
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Next,  the  poweron  flag  is  set  to  1,  and  InitxramO  is  called.  The 
images  are  displayed  by  sending  the  proper  flags  to  the  grabbers.  Then 
a  1024  byte  array  of  data  is  used  by  the  function  InitxramO  to  copy 
data  from  disk  to  "Xram",  the  grabber  control  ram,  and  located  at 
machine  address  0xa0000  (the  eleventh  block  of  2**16  bytes — beginning 
at  decimal  655360.  Grabber  ram  runs  from  655360  to  720695,  each  board 
has  the  same  addresses,  but  is  selected  by  a  flag  fed  to  the  board  via 
a  conputer  output  port.  Similarly,  if  parameters  are  to  be  changed, 
they  also  occupy  logical  address  space  beginning  at  0xa0000,  but 
another  flag  is  used  to  tell  the  boards  that  the  reference  is  to  the 
data  in  Xram. 

An  image  is  displayed  from  each  camera  on  the  monitor  and  the  computer 
waits  while  operator  adjusts  the  irises  of  the  cameras.  The  desired 
iris  setting  provides  an  average  light  flux  across  all  the  pixels  of 
128,  half  the  range  of  the  converted  number.  The  gain  and  offset  of 
the  A3C’s  can  be  adjusted  so  that  for  the  same  iris  opening  a  uniform 
sky  (either  all  clear  or  all  overcast)  each  camera  produces  the  same 
average.  With  the  lenses  in  the  test  configuration,  fl&  is  suitable 
from  two  hours  after  sunrise  until  two  hours  before  sunset  if  the  sky 
is  no;  heavily  overcast.  As  darkness  approaches,  f  stops  have  to  be 
rathe"  rapidly  decreased  (opening  the  iris),  and  after  the  last  stop  is 
reached,  time  exposures,  using  the  inject-inhibit  option  are  required. 
Mi  see  L  laneous  funct  ions,  Tot  1  ight  () ,  SetirisO.  and  SetiiO  have  been 
written  to  be  included  with  any  system  equipped  with  auto-iris. 

When  the  irises  are  set  to  the  satisfaction  of  the  operator,  any  key 
entry  initiates  the  endless  loop,  "while  (poweron*l)  {...>.  MainO 
proceeds  through  the  functions  listed  above,  ending  with  reportO  and 
then  Returns  to  the  beginning  of  the  endless  loop  executing  all  the 
functions  again. 

A3.  2.  4  GrabimageOt 

Initially,  when  the  exposure  is  satisfactorily  set,  any  key  can  be 
struck,  or  thereafter,  upon  return  from  the  end  of  the  loop  GrabimageO 
executes.  The  time  is  determined  from  the  on-board  clock  by  a  function 
"btimerO"  which  reads  the  clock  and  converts  the  hours,  minutes  and 
seconds  to  a  number.  Time2  and  time3  are  computed  and  stored  and  the 
addrebs  in  grabber  ram  into  which  the  pictures  are  to  stored 
calculated. 

.Initially,  the  grabbed  images,  those  being  observed  to  set  the  irises, 
are  scored  in  grabber  ram  at  address  0  up  to  the  size  of  the  image. 
The  first  true  image  is  collected  in  one  frame  after  the  address  is 
changed  by  Grab (a, b, c),  the  second  true  image  is  collected  in  one  frame 
and  the  address  increased  again.  At  the  third  address  change,  the 
image  is  turned  on  and  stays  on,  being  overwritten  into  this  fourth 
block  until  the  next  observation,  at  which  the  first  true  image  of  the 
second  observation  is  stored  ,n  the  second  block,  and  so  oh.  The 
reasons  for  these  gyrations  have  to  do  with  the  interactions  of  the 
various  flags  and  the  desirability  of  being  able  to  see  the  image  in 
between  observat ions. 
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The  TV  image  coming  in  from  either  camera  is  essentially  the  same  as 
any  TV  transmission  protocol.  At  ths  beginning  of  the  frame,  the  pixel 
in  the1  upper  left  corner  of  the  chip  is  read,  amplified  and 
transmitted.  Then  the  next  pixel  in  the  horizontal  row  is  read,  and  so 
on  until  all  360  pixels  in  the  row  have  been  read  and  transmitted.  '  The 
number  360  is  one  of  the  parameters  fed  to  the  grabber  control  table  by 
Init  xram < ) . 

Then  a  delay  occurs,  called  the  horizontal  blanking  time,  while  the 
putative  TV  receiver  moves  its  beam  back  to  the  first  pixel  of  the  next 
row,  after  which  the  second  row  is  transmitted  in  the  same  way.  This 
process  continues  until  all  240  rows  have  been  transmitted,  at,  which 
time  a  longer  delay,  called  the  vertical  blanking  time  occurs,  during 
which  the  beam  is  returned  to  the  upper  left  corner  of  the  chip. 
During  transmission,  the  chip  has  been  gathering  light  for  the  next 
picture. 

After  conversion  to  a  number  in  the  computer's  8-bit  flash  ADC,  the 
pixel  values  are  stored  in  grabber  memory  cell  <x)  where  x« (row 
number) * (pixels  per  row)+  column  number.  x  is  actually  a  single  cell 
address  which  increments  by  one  with  each  new  entrant,  and  starts  at 
the  address  provided  to  the  grabber.  This  cell  can  be  anywhere  in 
grabber  ram,  but  its  value  must  be  transmitted  to  the  grabber,  and  time 
must  be  allowed  to  change  the  address.  The  horizontal  blanking  time  is 
not  long  enough  for  this,  but  the  vertical  is,  hence  the  process 
requires  that  the  beginning  of  vertical  blanking  be  identified,  the 
input  turned  off,  the  address  changed  during  this  interval,  and  then 
the  input  turned  on  again.  That  is,  the  storage  of  numbers  coming  in 
must  begin  with  the  first  number  after  vertical  blanking  is  finished 
and  terminate  for  this  block  of  memory  when  vertical  starts.  Then  the 
address  can  be  changed  and  the  next  block  of  data  stored  somewhere 
else.  There  is  room  in  the  ram  for  12  complete  86400  vaxue  blocks,  but 
they  must  not  be  allowed  to  overlap,  nor  can  acquisition  be  started  in 
the  middle  of  a  frame  reliably. 

Thus  "grabO"  performs  four  tests.  If,  when  first  entered,  vertical  is 
off,  (vertical  retrace  is  in  progress)  it  is  not  known  how  long  this 
has  been  going  on  and  hence  there  may  not  be  time  enough  left  to  insert 
the  address  at  which  the  next  picture  will  be  started.  Hence  the 
machine  first  makes  sure  the  video  is  being  recorded  in  the  old 
address,  waits  until  vertical  comes  back  on  and  a  picture  is  being 
input.  A  bit  in  an  output  port  of  the  grabber  is  zero  while  input  is 
occurring.  At  the  end  of  the  picture,  that  bit  goes  to  1,  and  the 
vertical  blanking  interval  has  begun.  The  address  of  the  new  block  of 
grabber  ram  into  which  the  next  picture  is  to  be  placed  is  copied  into 
the  appropriate  registers  from  the  formal  parameters  a, b, c  in 
" grab (a, b, c) "  and  the  machine  then  waits  until  this  vertical  retrace 
interval  is  over  (about  1.6  ms.).  When  the  bit  drops  to  0,  the  machine 
tests  the  bit  continuously  until  the  vertical  retrace  begins  again — the 
bit  goes  to  1 — and  at  that  time  turns  the  video  off,  so  that  no  more 
data  will  be  recorded  until  the  address  is  changed  to  the  next  block. 

The  address  of  the  bloc'c  into  which  the  frame  is  placed  increments  by 


Micro  Science  Inc. 


Page  £6 


AFOSR  F49620-86-C-0093 


86400  bytes  per  image,  beginning  at  86400.  This  is  hen  5460  times  4  (a 
2  bit-offset  is  required  in  the  address  field. )  Hence  the  addresses  are 
manufactured  as  a  triple  of  2-digit  hex  numbers,  aabbcc,  one  pair  to 
each  of  three  ports,  where  addr  starts  at  21600  and  increments  by  this 
amount  each  frame,  a=addr/65536,  d=addr-a*65536,  b=d/256,  c=d-b#256, 
all  unsigned  integers  or  long  (32-bit)  integers.  Thus  the  first  four 
block  addresses  are  0x5460,  0xabc0,  0xfd£0  and  0x15180  of  which  only 
the  first  three  are  currently  used,  and  the  hex  values  of  a,  b  and  c  are 
respectively  00,54,60,  then  00, ab, c0,  then  00, fd,£0  and  finally 
01,51,80.  The  parameter  nsamples  is  currently  3,  but  can  be  increased 
to  11  if  needed  to  allow  a  larger  number  of  pictures  to  be  taken  at  any 
given  observat ion.  For  complex  or  slow  moving  skies,  this  may  prove 
useful. 

Each  board  is  controlled  by  the  same  signals  from  the  computer,  but  the 
data  are  collected  independent ly  in  each  grabber  ram  from  its  own 
camera. 

After  the  first  sample  is  collected  by  "grabimage ( ) ",  the  interval 
between  samples,  normally  2  seconds,  is  added  to  the  time,  and  the 
clock  is  examined  continuously  until  that  time  is  reached.  Then 
another  sample  is  collected.  When  nsamples  have  been  collected, 
control  returns  to  MainO  which  calls  GratozimgO.  There  is  a  possible 
error  in  this  protocol,  but  it  appears  to  have  no  effect  oh  the  output. 

A3.  2.  5  GratozimgO: 

The  job  of  this  function  is  to  copy,  from  each  board,  2  samples  of 
86400  bytes  each  through  a  65536  byte  window,  into  computer  ram,  break 
up  the  resulting  blocks  into  four  21600-byte  long  segments  each, 
transform  the  data  items,  which  are  8  bit  numbers  (characters)  into 
integers  (15  bits  plus  sign)  and  copy  a  selected  subset  of  the  data 
into  2  tables,  zimg0  and  zimgl. 

The  problem  lies  in  the  fact  that  the  window  "slides"  along  the  grabber 
ram  memory  space  65536  bytes  at  a  time,  while  the  data  from  the  cameras 
has  been  blocked  into  grabber  ram  in  86400  byte  contiguous  chunks.  The 
first  chunk  begins  at  address  86400,  which  extends  86400-65536=20864 
bytes  beyond  the  end  of  the  first  window.  So  the  window  is  shifted  to 
65536,  and  everything  from  window  +  20864  to  window  +  20864  +  21600  -  1 
copied  to  a  block  of  computer  ram  called  img0Cl.  Then  the  block  from 
42464  to  64063  is  copied  into  the  same  area  and  transfered  to  zimg0  by 
the  function  CopyimageO. 

There  are  now  only  1472  bytes  remaining  that  can  be  reached  through 
this  window,  so  these  are  copied  into  the  low  end  of  img0.  Then  the 
window  address  is  incremented  by  21600  and  the  remaining  20129  bytes 
brought  into  the  high  portion  of  irnglCl  from  the  low  end  of  the  window, 
and  transferred,  again  by  CopyimageO.  This  process  is  repeated  for 
both  boards,  until  16  records  of  21600  bytes  each  have  been  input. 

A3.  2.  5.  1  CopyimageO 

This  function  transfers  only  the  data  in  quadrants  1  and  2  (the  middle 
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two)  and  only  the  data  from  times  t  =  1  and  t«2,  and  of  these,  only 
alternate  rows  and  and  only  the  alternate  columns  from  30  to  298,  the 
middle  150  pixels.  Thus  the  88400  byte  picture  is  reduced  to  9000 
integers  by  CopyimageO  for  each  time,  far  each  board.  The  first  valid 
picture  is  stored  in  zimg0Cboard] Crow! Gcolumnl  and  the  second  picture 
in  zimgl  C board  1  Crl  Ccolumnl .  The  cast-outj  data  does  decrease  precision, 
but  adjacent  pixels  appear  to  have  some  qort  of  electronic  bias  which 
produces  a  "comb"  effect.  This  may  be  a  left-over  bit  in  the  ADC. 
Further  experiments  will  be  required  to  determine  the  trade-off. 

A3.2.S.2  Zimg  (t )  out  ( )  : 

After  the  two  pictures  are  taken  and  transferred  to  zimg Ct3 Cb] Cr3 Cel,  2 
36000  integer  arrays  are  in  memory.  These  are  filed  as  4  blocks  of 
18000  integers  to  a  harddisk  file  rcp\zimg2. 

A3.  2.  5.  3  ShowcloudO: 

Next,  if  desired,  the  service  routine  ShowcloudsO  executes,  retrieving 
the  files  just  stored  transforming  the  integers  to  characters^  and 
re-filing  them  into  4  21600-byte  arrays  Which  are  concatenated  into  a 
new  file  RCPZIMG3. 000,  and  returned  td  hard  disk,  in  which  form  they 
can  be  retrieved  by  the  SVPCIP  "file  read"  instruction  and  displayed  on 
the  monitor,  providing  operator  with  the  image  of  all  quadrants  that 
are  subsequently  to  be  processed. 

A3. 2. 6  Zolay0(zlmg0, zofs0) ,  Zolayl < zimgl , zofsl ) * 

At  this  point,  all  the  pre-processirig  of  the  camera  image  has  been 
completed,  and  all  the  data  to  be  processed  for  this  observation  are  in 
two  arrays  of  integers,  zimg0Cboardl Crowi [column!  and  and  zimgl C3 Cl CD, 
for  the  two  times.  First  the  layers  are  found  by  overlaying 
zimg0C0inC3  onto  zimy0Cl J Cl  Cl ,  by  routine  Zolay0  for  time  0,  and  the 
same  layers  found  again  at  time  1  using  zimgl C01 Cl  Cl  and  zimgl Cll Cl  Cl. 
The  two  routines  are  identical  except  fori  the  arrays  involved. 

j 

First,  all  the  values  of  zimgtC01  and  zilmgtCll  arc  summed  into  two 
numbers  meena  and  meenb  where  a  and  be  refer  to  the  boards.  Then  the 
normalizing  number  dmeen  equal  to  meena/meenb  is  made.  Every  cell  in 
board  b  is  multiplied  by  dmeen,  resulting  in  two  data  sets  in  zimg0C01 
and  zimg0Cll  whose  means  are  identical.  The  purpose  of  this  exercize 
is  to  compensate  for  the  irreducible  differences  in  the  two  ADC’s  or 
camera  voltages  which  do  vary  over  time  Sven  if  the  gains  and  offsets 
of  the  two  grabbers  are  tuned  ver^y  closely.  The  reason  for 
normal izat ion  is  to  assure  that  if  two  pixels  are  of  the  same  point  in 
the  sky,  the  only  difference  between  them  will  be  related  to  parallax 
effects,  not  by  camera  effects.  Since  i  the  former  appears  not  be 
significant,  the  difference  is  often  very  close  to  zero,  enhancing  the 
accuracy  of  the  overlay  process. 

The  next  step  is  the  actual  overlay  process.  Since  for  triangulat ion, 
negative  overlay  is  meaningless,  zimgtC01  is  placed  directly  over 
zimgtCll  and  shifted  "right"  one  column  at  a  time,  the  number  of 
columns  shifted  being  the  "offset".  As  the  shift  proceeds,  the  sum  of 
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the  squares  of  the  differences  between  corresponding  pixels  in  the  two 
images  will  change,  depending  on  the  number  of  pixels  in  the  overlay 
that  correspond,  that  is,  report  the  light  received  from  the  same  point 
in  the  sky.  For  any  given  offset,  the  square  root  of  the  the  resulting 
number  divided  by  the  number  of  points  compared,  is  retained  in 
zofst Coffset 3 ,  and  the  minima  in  it  found  later.  These  minima  are  the 
layers. 

A3.2.7  Calcz  O  ; 

After  the  two  ZolayO's  have  run,  two  lists  of  150  numbers,  zofs0C3  and 
zofslCl  are  in  memory.  The  next  step  is  to  find  the  minima  and  clean 
up  the  lists  of  minima  so  that  as  nearly  as  possible  the  layers  found 
at  t«0  are  matched  with  the  layers  found  at  t*l.  The  first  step  is  to 
smooth  the  zofs  themselves  so  the  minimum  finding  algorithm  works 
correctly.  This  has  been  done  in  two  ways}  either  the  zofs  list  can  be 
smoothed  with  a  Fast  Fourier  smoothing  routine,  or  each  point  can  be 
set  to  the  average  of  the  five  adjacent  points  of  which  it  is  the 
central  value.  Both  schemes  work  well.  The  former  is  more  accurate 
(returns  a  better  estimate  of  the  true  value  of  the  point,  assuming 
that  the  curve  is  smooth),  while  the  second  is  somewhat  faster,  and 
takes  less  program. 

A  minimum  is  defined  as  any  point  in  the  set  of  zofs  which  is  equal  to 
or  less  than  its  predecessor  and  less  than  its  successor.  If  the  curve 
has  a  U-shaped  bottom  which  could  occur  due  to  noise,  it  is  possible 
for  "minima"  to  be  found  two  points  apart  by  this  rule.  The  filter 
technique  essentially  eliminates  W* s  and  thus  eliminates  these 
occasional  strings.  On  the  other  hand,  the  cleanup  routine  matches 
points  which  are  most  nearly  equal,  so  that  if  two  nearby  points  are 
found  by  the  minimum  finding  routine,  for  both  curves,  two  minima  will 
remain  only  if  both  sets  produce  the  same  pairs,  somewhat  unlikely. 
Since  the  zofs  are  double  length  floating  point  numbers  the  equality 
test  is  not  necessary. 

Uhen  a  minimum  is  found,  its  index  is  placed  in  array  dwCUCt],  the 
value  of  zofs  at  this  point  placed  in  array  zdwCllCtl  and  1,  the  layer 
number,  is  incremented.  If  1  reaches  the  maximum  size  of  the  array, 
the  largest  minimum  in  the  array  is  deleted,  and  all  the  minima  below 
it  in  the  list  pushed  up  one.  In  a  reasonably  well  behaved  sky,  this 
is  unlikely  to  happen,  since  the  size  of  the  array  is  three  times  as 
large  as  any  reasonable  number  of  layers  (3).  On  the  other  hand,  it  is 
not  yet  known  what  the  consequences  will  be  of  a  towering  cumulus  whose 
perpendicular  face  is  visible  to  both  cameras.  That  could  produce  a 
lot  of  layers,  so  that  exceeding  array  size  is  an  indicator  of  a 
special  situation,  or  cloud  type,  no  necessarily  of  signal  noise. 

After  both  strings  are  processed  and  the  minima  found,  the  stacks  are 
processed.  A  square  array  is  filled  in  with  the  absolute  values  of  the 
differences  of  each  layer  value  in  one  board  with  each  layer  value  in 
the  other.  Then  any  row  colunmn  pair  intersection  which  contains  a 
zero  is  flagged  in  both  row  and  column.  Then  each  intersection  which 
contains  a  1  and  is  not  flagged  in,  either  row  or  column  is  flagged  in 
both.  Then  each  row  column  intersection  which  contains  a  2  and  is  not 
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flagged  in  either  row  or  column  is  flagged.  This  process  reduces  the 
lists  to  two  sets  of  equal  length  which  contain  at  equal  indices  layers 
of  equal  or  nearly  equal  offsets. 

Offsets  in  one  list  which  have  no  match  within  2  in  the  other  are 
discarded.  Each  list  is,  of  course,  the  layers  found  by  triangulat ion 
from  both  cameras  at  one  time.  The  reason  for  this  clean  up  is  that  in 
the  next  step,  the  layers  must  be  taken  in  pairs  of  nearly  equal 
offsets  for  the  two  different  times  in  order  effectively  to  measure 
velocity. 

A3. 2. 8  Oxdyl < ) t 

Beginning  with  the  lowest  layer  (the  largest  offset),  the  images  from 
board  0  for  times  0  and  1  are  overlaid  for  correspondence.  If  for 
both,  zdw  (the  value  of  zofs  for  this  layer)  is  zero  (no  layer  has  been 
entered  in  this  slot)  or  the  value  of  the  offset  is  greater  than  140 
(the  layer  is  too  low  to  be  have  a  suffiently  large  number  of  points  to 
be  reliable)  the  layer  index  is  decremented  and  the  next  layer 
examined.  It  is  possible  for  there  to  be  no  layers  at  all,  in  which 
case  the  words  "clear"  are  printed  in  the  report. 

When  a  layer  is  found,  the  first  step  is  to  make  an  image  restricted  to 
just  the  points  in  this  layer,  so  that  just  these  can  be  overlaid. 
Function  Makevimg (ofs0, ofsl, layer  number, zimg0, zimgl,  vimg,  cover)  is 
used  for  this  purpose.  (The  string  of  labels  following  the  function 
name  are  the  formal  parameters  which  tell  the  function  which  is 
compiled  and  in  the  library,  the  value  of  the  two  offsets  (one  from 
each  list  for  the  same  layer,  very  likely  equal),  the  layer  number  and 
the  locations  of  the  four  arrays  containing  the  data  Makevimg  will  use 
or  make. ) 

A3. 2. 8.  1  Makevimg () i 

The  function  treats  the  two  zimgs  in  turn.  First  a  normalization 
number,  dmean,  is  found  for  zimg0Cll C i 3 Cj+ofs01  in  terms  of 
zimg0C01 til C j 3  for  the  area  which  overlaps  (which  may  be  different  from 
the  total  image).  Then  a  zero-level  for  this  layer  is  found  as  the 
mean  of  zimg0  for  the  area.  This  value  is  copied  to  all  the  cells  in  a 
new  array  vimgCtl Cil Cjl  which  are  not  in  the  overlapping  area.  Then 
the  standard  deviation  of  all  the  points  in  the  layer  is  found  as  the 
square  root  of  the  average  of  the  sum  of  the  squares  of  the  deltas  (the 
difference  between  the  values  at  one  time  and  the  corresponding  values 
at  the  second  time). 

Then  the  deltas  are  computed  again  (they  should  have  been  saved,  but 
space  is  tight),  and  if  a  delta  is  within  one  standard  deviation  of  the 
mean  of  the  array  the  value  of  that  pixel  from  zimgt  C0"1  Cil  Cjl  is  copied 
to  an  array  vimgCtl Cil Cjl,  and  (for  t*0  only),  a  counter,  coverClayer 
number],  is  incremented.  Otherwise  the  zero-level  is  copied  to  that 
cell  and  cover  is  not  incremented.  Cover,  therefore  at  the  end  of  this 
exercise  is  a  count  of  all  the  points  in  this  layer  that  correspond  and 
when  divided  by  the  number  of  points  in  the  area  considered  is  in  some 
sense  a  measure  of  the  share  of  total  cover  provided  by  this  layer 
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alone.  The  pixel  value  should  be  compared  with  a  cutoff  value, 
representing  blue  sky,  and  the  point  excluded  if  it  falls  below  this 
cutoff.  This  is  not  done  at  present,  hence  cover  over-represents  the 
actual  degree  of  cloudiness  by  not  picking  up  ''holes"  in  the  overcast 
if  they  exist. 

The  same  process  is  performed  on  zimgl C03 Cil Cj, ofsll  to  produce 
vimgCll Cil Cjl  and  two  arrays,  of  the  same  layer  of  the  overlapping  part 
of  the  same  sky  at  two  different  times  are  ready  for  analysis.  DxdylO 
proceeds  to  move  one  image  around  the  other,  first  with  the  center  of 
one  image  exactly  on  the  center  of  the  other,  using  the  smaller  of  the 
widths  of  the  images  as  the  limiting  width  of  the  two.  Again,  using 
this  time  a  function  DvolayO,  the  deltas  are  computed  for  the  points 
in  the  common  area,  summed  and  averaged  to  produce  a  value,  vofs,  for 
this  particular  amount  of  shift  of  one  image  over  the  other. 

A3.  2.  8.  2  DvolayO: 

In  the  first  call  to  DvolayO,  there  is  no  shift  du=0, dv*0;  in  the 
successors,  one  image  is  moved  first  one  pixel  right  du=l,dv=0,  and 
DvolayO  is  called  to  find  vofs  for  the  overlapping  area  for  this 
shift.  Then  vimgC01  is  shifted  one  cell  up  (du  is  now  1  and  dv=l)  and 
vofs  calculated.  Then  the  shift  is  two  left,  then  two  down,  then  three 
right  and  so  on.  At  each  shift, the  vofs  produced  is  compared  with  its 
predecessor.  If  it  is  less,  it  replaces  the  minimum  and  the  shift  in  x 
and  y  for  this  new  minumum  vofs  saved.  If  after  the  last  minimum  was 
found  (successive  values  of  vofs  were  all  higher)  and  the  current  shift 
is  at  least  two  steps  outside  that  minimum  all  the  way  around,  then 
that  shift  is  accepted  as  the  true  local  minimum.  The  corresponding  du 
and  dv  of  that  minimum  are  now  directly  proportional  to  the  velocity  of 
this  layer  in  x  and  y,  that  is,  along  the  axis  of  the  cameras  and 
perpendicular  thereto. 

The  process  is  repeated  for  all  the  layers. 

A3. 2. 9  Report  < ) i 

Report  converts  the  values  of  dwClayer] , duClayerl, dvClayerl,  which  were 
found  in  the  preceding  steps,  to  actual  altitudes  and  directions,  and 
displays  them  on  the  monitor  for  the  user.  The  following  computations 
are  made. 


Altitude  in  feet*kl*k2/ (k3*dwCln3 )  +  k4  where 

kl=focal  length  of  the  cameras  in  pixel  widths  *720, 
k2*distance  between  the  two  cameras  =  963.52  feet, 
k3=number  of  pixels  per  point  in  the  image  *2, 
dwCln3=z-of fset  of  this  layer. 

k4  »  the  altitude  ASL  of  the  station  *  240  feet. 

Thus  altitude  of  the  layer  is  346867. 2/dwCln3  +  240.0  feet  ASL,  +/'-  2 

inches.  The  altitude  error  is  actually  at  least  +/-  1  pixel  width  at 
that  altitude.  If  dwCln3=94,  say,  then  the  reported  cloud  layer  would 
be  3690+240=3930.  But  94  is  the  integer  reported  for  any  dw  greater 
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than  93.5  and  less  than  94.5,  so  that  the  layer  lies  somewhere  between 
3910  and  3950  feet — one  percent.  More  to  the  point,  the  value  of  kl  is 
an  approximat ion,  arrived  at  as  follows: 

The  width  of  the  pixels  in  x  is  22.2  micrometers.  The  focal  length  of 
the  lens  is  IS  mm.  Then  kl*16/.  0222-720.  But  focal  length  changes 
very  slightly  with  iris  opening  and  certainly  with  focus  (the  cameras 
are  maintained  at  "infinity"),  the  pixel  width  is  an  engineering 
specif ication  which  undoubtedly  has  a  small  unknown  error,  and  the 
cameras  lines  of  sight  are  almost  certainly  not  parallel  nor 
perpendicular  within  some  small  error,  for  example  the  fact  that  the 
lines  to  the  center  of  the  earth  from  each  are  not  parallel. 

These  errors  accumulate  randomly  and  could  average  out  to  zero,  but  to 
be  safe,  in  the  worst  case,  we  suggest  that  measured  altitude  is  true 
altitude  within  5X.  If  the  cloud  base  in  the  range  of  view  is  not 
level,  the  altitude  produced  is  the  average  altitude,  about  that  over 
the  midpoint  between  the  two  cameras. 

After  altitude  is  found  for  each  layer,  velocity  and  direction  are 
computed  as  follows: 

dx*du*alt*k4/k5/k6*k7  where 

du  is  the  offset  returned  by  DxdylO  for  this  layer, 

alt  is  the  altitude  of  this  layer  in  feet  as  found  above 

k4  is  the  number  of  pixels  per  unit  of  x  offset  -2 

k5  is  focal  length  in  x-pixel  widths  »  720 

k6  is  the  interval  between  the  exposures,  2  seconds 

k7  converts  feet  per  second  to  knots*3600/6080 

k4*k5/k6*k7». 0008223687,  the  constant  found  in  the  code 

Next,  dy  is  found  in  the  same  way,  but  with  an  adjustment  for  the  fact 
that  the  y-dimension  of  a  pixel  is  1.17  times  the  x-dlmension.  .  This 
implies  that  unit  velocity  in  y  will  cover  less  pixels  than  the  same 
velocity  in  x,  so  that  k4(y)*  2*1.17*2.34.  Thus  the  constant  found  in 
the  code  is  .000962775.  The  other  coefficients  are  the  same. 

After  dx  and  dy  are  found  for  a  layer,  Report ()  computes 

Velocity  ■  sqrt (dx  2+dy  2),  knots. 

First,  if  cloud  movement  is  along  the  axes: 

If  dx*0  and  dy*0  then  direction  is  zero  (and  velocity  is  zero  as 
well). 

If  dx*0  and  dy) 0  then  direction  is  baseline  direction. 

If  dx*0  and  dy(0  then  direction  is  180  plus  baseline  direction. 

If  dy*0  and  dx)0  then  direction  is  90  plus  baseline  direction. 

If  dy»0  and  dx  <0  than  direction  is  270  plus  baseline  direction. 

Otherwise,  angle  *  arctan(dy/dx)*degrees/radian,  in  degrees  relative  to 

the  base  line  direction,  27.3037  true. 
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If  dx <0  than  direction  »  270+basel ine  -  angle, 

If  dx>0  then  direction  ■  90  +baselin*  -  angle. 

After  a  header  is  printed  giviny  the  name,  latitude  and  altitude  of  the 

station,  the  date  and  the  time,  column  headings  are  printed  as  "cloud 

base  ASL  moving  toward  deg  true,  velocity  kts,  percent  cover",  and  the 
altitude,  direction,  altitude  and  cover  are  printed,  one  row  for  each 
layer,  from  lowest  to  highest. 

"Report ()"  returns  to  main,  which  spins  until  time  for  the  next 
observation,  at  which  point  the  entire  program  repeats. 


A3. 3  Annotated  Code.  Attached  after  this  Annex 

The  operating  system  and  its  functions  which  run  in  interpretive  mode 
(RCP)  is  given  first,  followed  by  listings  of  FEBLIB. C,  the  compiled 
code,  called  from  RCP.  After  this,  listings  of  the  ancillary  programs 
which  were  developed  but  not  included  in  the  operating  code  are  given. 
No  proprietary  software  (to  my  knowledge)  is  included  except  the  WHSKY 
program  itself,  which  is  copyrighted  by  Micro  Science,  Inc. 

A3. 4  Miscellaneous  software 

Several  programs  not  included  in  the  operational  WHSKY  were  written 
during  the  development  process.  These  are  discussed  below  and  the  code 
is  attached  following  the  operational  system  code  (see  above). 

A3. 4.  1  Fourier  noise  filter,  FilterOi 

This  algorithm  was  devised  by  Aubanel  and  reported  in  Byte,  2/B4.  It 
has  been  translated  to  C  and  compiled,  and  appears  to  work  exactly  as 
specified.  It  is  interesting  to  note  that  the  routine  returns  the 
coefficients  of  the  frequencies  which  are  summed  to  produce  the  cleaned 
up  function.  If  a  time  series  of  cloud  pictures  is  made,  the 
coefficient  of  the  fundamental  frequency  displays  a  shift,  and  the 
frequency  of  this  shift  relative  to  the  fundamental  is  in  fact 
proportional  to  the  velocity  of  the  cloud. 

A3. 4. 2  Lens  compensat ion,  Adjust  lens () i 

If  the  lens  were  a  theoret ical ly  perfect  pinhole,  the  distance  between 
any  two  points  in  the  sky  at  the  same  altitude  would  be  exactly 
proportional  to  the  distance  between  them  on  the  plane  of  the  chip.  If 
the  points  were  projected  onto  a  sphere  rather  than  a  plane,  points 
away  from  the  perpendicular  would  appear  closer  to  the  center.  That 
is,  if  a  picture  of  a  perfect  grid  on  a  plan*  normal  to  the  line  of 
sight  were  taken  with  a  standard  wide  angle  lense,  the  grid  lines  at 
the  periphery  of  the  picture  would  appear  closer  together  than  those  in 
the  center,  resulting  in  the  "pin  cushion"  effect.  Thus  the  camera 
"lies".  The  error  is  small,  less  than  a  pixel  width  for  most  of  the 
field  of  a  16  mm  lens,  and  because  other  errors  are  far  greater,  the 
compensatory  routine  for  this  error  is  available  but  not  part  of  WHSKY. 
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A3. 4. 3  File  processing  routines! 

In  the  development  process,  it  is  highly  desirable  to  save  intermediate 
results  so  that  checked  out  routines  do  not  have  to  be  re-run  to 
produce  the  input  for  a  function  which  is  under  development.  Hence  it 
is  common  practice  to  store  these  intermediate  results  along  the  way. 
In  the  final  version,  which  need  save  nothing,  the  routines  have  been 
deleted. 
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Annex  4.  Source  code  listings 


/*  RCPNWHSKY. 800  */ 

/*  Complete  operating  system  8/18/87  */ 

/*  Copyright  (C>  1987  by  Micro  Science  Inc.  ALL  RI6HTS  RESERVED*/ 

/*  Development  supported  in  part  by  US  Air  Force  OSR  */ 

/*  Contract  F49620-86C-0093DEF  */ 

♦define  npin  21800 
#define  nl  8 
♦define  nsamples  3 

/•  npin“21600*60*360  one  fourth  of  a  grabbed  image  */ 

/*  nt  *  number  of  grabbed  images  per  observation*/ 

/*  nsamples,  one  more  than  required  for  dx,dy  at  a  given  obs’ n. */ 

/*  for  a  single  cloud  image  (8  files),  set  nsampless2  */ 

/*  arrays  should  be  redimensioned  to  C23C1203C1353  */ 

static  float  t ime, delayt ime; 

static  float  t ime2, t ime3, t ime4, dt 2=600, dt3=2, dt4*1200; 

/♦next  observation,  next  sample,  intensity  wait*/  , 

/*dt3  should  be  increased  if  altitude  of  layer  is  high  */ 
static  char  img0Cnpirt3 ; /*holds  grabber  part  for  copy image*/ 
static  int  zimg0C23 C603 C1503, zimgl C23 C603 C1503 ;  /*t,  brd,row,col  */ 
static  int  vimg C23 C603 C1503 ;  /*t ime, r, c <board0)  for  one  layer  */ 
static  double  zofs0C1503, zOfsl C1503,  vofsC3023, coverCnll ; 

/*z  offsets,  at  times  0, l,vofs  at  any  time,*  cover  */ 
static  int  dwCnl3 C23, zdwCnl3 C83,duCnl3,dvCnl3 ;/*by  lr,alt,vx,vy  */ 
report , */ 

mainO  < 

runcC'SHELL  SVPCIP  -t  0xa0000")f  /*initialize  grabber  */ 

rune < "LIBL  FEBLIB") ) 

runcCLIBL  RCP\MATHLIB\MATHLIB")  ; 

rune ( "LIBL  RCP\GRAPHLIB\GRAPHLIB" ) ; 

runcCSET  BREAK  SCREEN  (2,  0,  0,  0)  {SCREEN  (0,  0,  0,  0)  ;  " )  ) 
rune ("SET  END  SCREEN (2, 0, 0, 0) ; SCREEN (0,  0,  0,  0) ; 
int  poweron=*l,  i=0; 

initxramO}  /*  read  in  XRAM.  001,  copy  to  0xa0000,  both  boards,  */ 
outp (0x301, 0xac) ;  /*video  on  both  boards*/ 

printf ("Adjust  irises  in  cameras.  When  ready,  hit  any  key\n"); 
get char ( ) ; 

while  (poweron  »  1)  <  /*  forever  */ 

grabimageO;  /*grab  both  images  for  t0  and  tl*/ 

gratoz img ( > ;  /*  brings  images  into  c-ram,  re-pots  into  zimg0, 1*/ 
zimgoutO;  /*  2  arrays  for  showclouds — adj’d  values  of  mid-sky  */ 
showclouds ( img0, zimg0, zimg 1 ) }  /*  makes  file  to  show  all  four 
zimgs  on  monitor  */ 

/*  printf  ("end  of  showclouds.  Key\n")j  getcharO;  */ 

/*  zimginO;  /*reeovers  both  arrays  */ 

/*  printzimgOj  */ 

zolay0(zimg0,  zofs0) ;  /*overlays  zimg0C03  onto  zimg0C13  */ 

/*  zof 30out < ) ;  */  /*stores  zofs0  produced  by  zolay0  */ 
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zolayl  (zimgl,  zofsl ) ;  /*overlays  zimglC01  onto  ziMglCl!!  */ 

/*  zofsloutO;  */  /*stores  zofsl  produced  by  zolayl  */ 

/*  zofs0in()f  */  /*recovers  zofs0  for  calcz  */ 

/*  zofslinO;  */  /*recovers  zofsl  for  calcz  */ 

/*  printzofsO;  */ 

/*  screenzof s O ;  */ 

calcz  O;  /*makes  dw, zdw,  the  layers,  from  the  zofs  */ 
dxdylO;  /*  makes  vimgs  by  layer  and  zimgs,  and  du,dv  */ 

/*  printdws  < ) ;  */ ' 
report ( ) ; 

delaytime  *  time£;  spinO; 

>  /*  while  */ 

>  /*end  of  whsky  */ 

/•  **##*#*##****#**#*###*  subroutines  ####*####**###*####*#***#*#**  */ 

btimerO  <  /*produces  the  time  in  seconds,  floating  */ 
char  acC103;  char  *acpoint;  int  ret}  float  hr, min, sec; 
dostime(ac, 4) 5  retasscanf (ac, "Xf :Xf :Xf", &hr, Amin,  tsec)  | 
time  «  3600*hr  +  b0*min+sec } 

/*  printfC'Time  in  btimeris  X.  2f \n",  t  ime)  }  getcharO;*/ 
ret  urn  < ) ;  > 

spinO  <  do  btimerO  |  while  (time  <  delaytime)}  > 
initxramO  < 

int  i, blent}  long  a»0xa0000;  char  xramCnuml}  char  *destinwxram} 

FIL.E  *inl ;  if  ( <  inl»fopen ( "RCPNXRAM.  001",  "rb")  >  l-NULL) 

<blcnt»fread (dest in,  num, 1, ini ) } f close  (inl);> 

else  <printf  < "f i leopen  RCPXXRAM. 001  f ai led" >  (exit  O  }> 
outp (0x301, 0x24) 1  /*ready  boards  (at  source)  to  receive  prams  */ 
runcO'SET  TRUST  ON")}  char  *source* (char  #)a| 
movmem  (dest in, source, num) ;  /*copy  prams  to  boards*/ 
runcO'SET  TRUST  OFF")}  o</tp(0x301, 0xac)  }  /*return  brds  to  digitg*/ 

> 

grabimageO  -C  /*moves  image  memory  into  SVPCIP  ram  for  gratozimg ( ) */ 
int  i«=0,  ck»-2,  v=*2;  long  addr, t 16*65536L ;  unsigned  a,  b,  c,  d; 
outp (0x301, 0xac> ;  /*camevas  on,  both  boards  */ 

btimerO;  timec^time;  /*  converts  dostime  to  a.  single  floater  */ 
t ime3“t ime2f  time2*time2  +  dt2;  addr«npin}a«®} b=0x54}c*0x60;d«0} 
for  ( i*=0; i (nsamples ; i++)  <  /*nsample*  from  cameras  */ 
print  f  (  "adr»Xld  i=Xd,  a**Xx,  b=Xx ,  c“Xx\n",  addr,  i,  a,  b,  c)  ; 
grab (a, b, c) ;  /*from  FEBLIB.  Finds  vert  blanking  time  (1.6ms)*/ 
addr+=npin;  a™addr  /  tl6;  /*calc  new  xram  memory  address  */ 
d«addr  -  a*tl6;  b“d  /  256}  c=d  -  b  *  256} 

time3+«dt3;  delayt  ime*t  ime3;  spin( )  }  >  >  /*  end  of  grabimageO  */ 
gratozimg  O  < 

/*  for  each  t,  for  each  grabber,  copies  4  panels  of  npin**£1600 
tytes  from  grabber  to  ram,  then  extracts  middle  alt  bytes 
integerizes,  copies  to  zimg (t ) CbD C60D [1503,  replacing 
copyimage () +makezimg ( ) .  zimgout  stores  to  disk. 

Cloud  can  be  displayed  on  the  monitor  under  SVPCIP  using  the 
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"showcloud < ) "  result.  Perform  f  r  RCPZIMG3. 000.  Note  no  \.  */ 

FILE  *out; 

char  nameC20D;int  a,  b,  t,  i,  im,  .1,  jm,  blent,  board,  quad,  ww*0; 
unsigned  delt=0; 

char  *source,  *destin*4:..g0,  *pl £73, *p4C6D  ; 

•pl^rcpNimg"  ;*p4«"002.  00" ; 

long  xram=0, window=0xa0000, tl£=65536L, upper»window+t 16; 
runcC'SET  TRUST  ON");  out  p  (0x308,  ww)  ; 
for  (board=0;board <2;board++>  \ 

if(board*=0)  outp(0x301, 0x28) ;  else  outp(0x301, 0x0c) ; 

nindow=0xa0000 ; ww=0 ; 

for  (t*0;t (nsamples;t++)  < 

printfCt  q  b  window  ww  delt  source  destin  sink  xramNn"); 
for  (quad*0;quad (4;quad*+)  { 

printf("X2d  Xd  Xd  X71d  ",  t, quad,  board,  window)  ; 
if  (.window  ♦  21600L  >  upper  >  < 

delt=*upper  -  window; source^ (char  *) window; 
printfC’  X2d  X6u  X71d  X71d  X71dXn", 

ww, delt, ( long) source, ( long) dest in, xram) ; 
movmem (source,  destin,  delt);  /*wove  the  bottom  part  */ 
ww++ ; out p (0x306, ww) ;  xram+*npin;  /amove  the  window  up*/ 
source*(char  *)0xa0000;  dejtin-(char  *>(img0  +  delt); 
printfC  X2d  X6u  X71d  XT^d  X71d\n", 

ww,  npin-delt, ( long) source, (long)destin, xram) ; 
movmem (source, dest in, 21600-delt ) ;  /*move  the  top  part*/ 
window*0xa0000  +  npin  -  delt;  /*finish  moving  the  window*/ 
dest in*img0;  > 

else  {source* (char  *)window;  movmem (source, dest in, npin) ; 
window*window  +  21600L;  xram+»21600;  printfC  X2d  ",ww); 
printfC  X71d  X71d  X71dXn", 

( long) source, (long) dest in, xram) ; > 

/•  one  block  of  char  raw  data  is  now  in  dest in**img0C2 16003. 

T*  q  is  0  or  3  or  t*0,  continue  (the  data  is  not  of  interest), 
else  write  it  to  a  file,  copy  alternate  bytes  to  zimg0C23 £603 Cl 503, 
zimgl £23 £603 Cl 501 — t,b,r,c.  Then  using  zimgouts  file  the  2  zimgs  */ 

if  (quad!^0  &&  quad!=3  &&  t>0)  { 

copyzimg ( img0, zimg0, zimgl, t-1, quad, board) ; 

>  /*if*/ 

>  /*quad*/ 

>  /*t*/ 

>  /*board*/ 

outp (0x301, 0xac> ;  /*video  on  both  boards*/ 
rune ("SET  TRUST  OFF"); 

/*  zimg0  and  zimgl  are  now  full,  ready  to  be  stored  by  zimgout*/  > 
/*  end  of  yratoimgO  */ 


calczO  {  /*from  zofs0,  zofsl,  finds  dw£ln3  £03,  dwClnl  £13  */ 
int  i, j, ka£23, k, 1, dt Cnl+13 £nl+13 ;  /*  first  5-point  smooth*/ 
vofs  £23  =  ( zof S0C03  +  zofs0  £13 +zofs0  £23  +zof s0£33+zof 50  £43 ) *0. 2L ; 
for  <j=3;j - 148; j++) 
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vofsCj3»  <5*vofsCj-13-zofs0Cj-33+zofs0Cj+23 ) *0. 2L; 
for  <j«2; J <148;j++>  zofs0Cj3=vofsCj3 ; 

/*  printzofsO  ;*/ 

zofs0C13**zofs0C23-l  ;zofs0C03*zofs0C13-l ; 

vofsC23“< zofsl  C03+zofsl  C13+zofsl  C23+zofsl  C33+zofsl  C43)*0. 2Lj 
for  (j«3;j <148|j++) 

vofsCj3-<5*vofsCj-13-zofslCj-33+zofslCj+23)*0.2L; 
for  <J=2; j <148; j++)  zofsl Cj3=vofsCj3 ; 
zofsl fl3*zof*l C23-1 ; zofsl C03=zofsl C 1 3  —  1 ; 

/*  screenzofsO  ;*/ 

for-  <j«l,  kaC03«0,  kaC13=0;j  <146;  j++)  < 

if  <zofs0Cj-13>  a*zofs0C  j3  44  zof s0Cj+l 3 > zof s0Cj3 ) 

{dwtka C03 3  C03 ; zdwCkaC03 3  C03 ■ < int ) zofs0Cj3 ;  kaC03++; 
if  <kaC03 > *nl )  <pushw<0> ;kaC0D*nl-l ; >  > 
if  <zofsl Cj-13> «zofsl Cj3  44  zofsl Cj+13> zofsl Cj3 > 

<  dwCkaC13 3  C13*j  ; zdwCkaC13  3  C13*  ( int )  zofsl  Cj3  ;  kad3*+j 
if  <kaC13>«nl)  {pushw < 1 > ; kaC13=nl-l ; >  >  > 

/*  clean  up  dw/zdw  stacks  */ 

for  <j*0; j <nl ; j++)  printf <"zdw, dws  %5d  *5d  X5d  X5d\n% 

dwCj3 C03, dwCj3  C13, zdwCj3 C03, zdwCj] C13 ) ; 
for  <i*0;  i  <nl ;  i++)  {  dt  Ci3  Cnl3*0;dt  Cnl3  Ci3*0;  /*initialize  flags*/ 
for  < j“0; j <nl ; j++)  {dt Ci3 Cj3«dwCi3 C03-dwCj3 C13 ; 
if  <dtCi3  Cj3  <0)  dtCi3t.)3  —  dtCi3Cj3; 

printf  <"*3dM,dtCi3Cj3)  ;  >  printf  <M\n“>  ;  >  /*fill  it  in*/ 
for  <k“0;k<3;k++)  {  for  ( i=0;  i  <nl ;  i++)  {  for  < j*0; j  <nl ;  j-*-*)  { 
if  <dt Cnl3 Cj3  ! »0  II  dt Ci3 Cnl3 !*0)  continue; 
if  <dtCi3Cj3--k)  <dt Ci3 Cnl3»l ;  dt Cnl3 Cj3-1 ; >  >  3  > 
for  <i«0, k»0; i <nl ; i++)  {  /*  push  up  ths  selected  i tarns  */ 

if  <dt  C i 3 Cnl3  >  0)  {dwCk3 C03»dwCi3 C03 ;zdwCk3 C03-zdwCi3 C03 ;k++; >  > 
for  <i-k;i <nl;i++>  {dwCi3 C03-0;zdwCi3 C03-0; >  /*  end  fill*/ 
for  <i-0,  k-0;  i  <nl ;  i++)  {  /*  push  up  the  sal  act  ad  items  */ 

if  <dt  Cnl3 Ci3>0)  <dwCk3 C13=dwCi3 E13 ;zdwCk3 C13-zdwCi3 C13 ;k+*;>  > 
for  <  i*k;  i  <nl ;  i-*-+)  <dwCi3  C13”0;zdwCi3  C13*0;  >  /*  and  fill*/ 
for  < j"0;  j  <nl ; j++)  printf<M  *5d  %5d  *5d  %5d  keyW, 
dwCj3  C03,dwCj3  C13, zdwCj3 C03, zdwCj3  Cl 3  > ; 

/•gatchar < ) ;*/  > 

pushw<b)  int  b;  {  int  i , imx=0,  mx-zdwC03 Cb3 ; 

for  < i»l ; i <nl ; i++)  {if  <zdwCi3 Cb3>mx)  {mx=zdwCi3 Cb3 ; imx»i ; >> 
if  <imx<nl-l)  <  for  < i»imx ; i <nl-l ; i++) 

<zdwCi3 Cb3«zdwCi+13 Cb3 ;  dwCi3 Cb3=dwCi+13 Cb3 ; >  >  > 


dxdy 1 ( )  /*v.  2/18/87  */ 

/*  dithering  routina — moves  center  of  vimg0  around  canter  of  vimgl*/ 
{ 

int  dul-0, dvl“0, mdu-0,  mdv»0,  t,  n,  1 ; 
int  ils, i2s, iend, jls, j2s, jend; 
int  ofs0, ofsl, jrng, f 1, f2; 
double  vof,minvofs; 

for  <  l-=nl-l ;  1>  =0;  1  — )  /*take  a  layer*/ 

if  <zdwCl]C03>0  44  dwC13C03  <140  44  zdwC13C13>0  44  dwmC13<140) 

{  ofs0=dwC13  C03 ;ofsl=dwC13  C13 ;f 1=0; f2=0;dul*0;dvl=0;mdu=0;mdv=0; 
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Makevimg (ofs0, ofsl, I, zimg0, zimgl, vimg,  cover) ; 

if  (ofs0>ofsl)  jrng»150-ofs0;  else  jrng*150-ofsl ; 
minvofs»dvolay (vimg,  0, 0, 60, 0, 0,  jrng) ; 

printf ("«invofs*Xlf  dul=*3d  dvl=S3d  jrng»*3d  cover®)  lf\n", 
minvofs, dul,  dvl, jrng, cover C 11 ) ; 
for  (t«l ;t (30;t+*2)  -C 

for  (n>=l  ;n  (=*t  ;n++)  <  dul++; 

if  (dvl)0)  < i ls*0; i2s=dvl ; iend=60-dvl ;  >  1 
else  { i ls»-dvl ; i2s=0; iend=60; > 
if  <dul>0)  <  j  ls=0;  j2s*=dul ;  jend=jrng-dul ;  >  ng;> 
else  <jls**-dul ;  j2s*0;  jend=jrng ;  > 
if ( (vof«dvolay (vimg, i Is, i2s, iend, j Is, j2s,  jend) ) (minvofs) 
<einvofs“vof ;mdu“dul ;mdv*dvl ; > 
printf  ("t«X3d  n«* 3d  dul«=X3d  dvl=X3d  vof»*lf  \nH, 
t, n, dul, dvl, vof ) ;  >  /*n*/ 
for  (n**l  ;n  («t  ;n++)  <  dvl  — ; 

if  (dvl)0)  < i ls=0; i2s=dvl ; iena»60-dvl ;  >  > 
else  -C  i  ls«*-dvl ;  i2s=0;  iend=60;  > 
if  (dul)0)  -C  j  ls=0;  j2s=dul ;  jend=jrng-dul ;  >  ng;> 
else  -Cjls=-dul  ;j2s=>0;  jend=jrng;  > 
if ( ( vof *dvolay (vimg, i Is. i2s, iend,  jls,  j2s,  jend) ) (minvofs) 
■Cminvofs«vof  ;mdu»dul  ;mdv=dvl ;  > 
printf  ("t*)t3d  n*X3d  dul=*<3d  dvl"*3d  vof**%lf\n  ", 
t, n, dul, dvl, vof ) ; >  /*  n  */ 
if  (mdv>dvl  +  l  &&  mdu(dul-l)  f  1=1 ;  else  f2*=0; 
printf ( "mdu*%3d  mdv*X3d  fl*X2d  f2«*%2d\n", mdu, mdv, f 1, f£) ; 
if  (fl— 1  &&  f2«-l)  break; 
for  (n*l ;n(*t+l ;n++)  {  dul — ; 

if  (dvl)0)  <ils«0;  i2s*»dvl ;  iend*60-dvl ;  >  > 
else  {i ls«-dvl ; i2s*0; iend=60; > 
if  (dul)0)  < j ls«0; j2s=dul ; jend*jrng-dul ; >  ng;> 
else  -Cj  ls*-dul ;  j2s*0;  jend=jrng;  > 
if ( (vof»dvolay (vimg, i Is, i2s, iend, jls,  j2s,  jend) ) (minvofs) 
•Cminvofs**vof  ;mdu=dul  ;mdv*dvl ;  > 
printf ( "t*X3d  n*=%3d  dul^Sd  dvl=*3d  vof*Xlf\n  ", 
t, n, dul,  dvl,  vof)  ; 

>  /*  n  */ 

for  (n*l ;n (=t+l ;n++)  <  dvl++; 

if  (dvl>0)  < i ls«0; i2s=dvl ; iend=60-dvl ;  >  > 
else  < i ls»-dvl ; i2s=0; iend=60; > 
if  (dul)0)  •(  j  ls=0;  j2s=dul ;  jendajrng-dul ;  >  ng;> 
else  {j ls»-dul ; j2s=0; jend=jrng ; > 
if ( ( vof =dvolay (vimg, i Is, i2s, iend, jls, j2s, jend) ) {minvofs) 
<minvofs=vof ;mdu=dul ;mdv=dvl ; > 
print f ( "t=*3d  n»K3d  dul=%3d  dvl=%3d  vof=%lf\n  ", 
t,  n, dul, dvl, vof) ; 

>  /*n*/ 

if  (mdv(dvl-l  mdu)dul+l)  f2*l;  else  fl=0; 

printf  ("mdu«=*3d  mdv=*3d  fl=*£d  f2=%2d\n",  mdu,  mdv,  f  1,  f2)  ; 
if  <fl«l  f2*=l )  break; 

>  /*t*/ 

d u C 1 1 =mdu ; dv C 1 ] =mdv ; 

>  /*  layer  */ 
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makevimg (ofs0, ofsl, 1, zimg0, zimgl, vimg,  cover) ; 

if  (ofs0)ofsl)  jrng»150-ofs0;  else  jrngs»150-ofsl ; 
roinvofs*dvolay (vimg, 0,  0, 80,  0, 0,  jrng) ; 

printf  ("minvofs*%lf  dul=*3d  dvl=*3d  jrng*X3d  cover*%lf\n", 
minvofs, dul, dvl, jrng, coverCll ) ; 
for  <t*l;t (30;t+*2)  < 

for  (n«*l  ;n  (=t  ;n++)  {  dul++; 

if  (dvl)0)  -C  i  ls=0;  i2s=dvl ;  iend=60-dvl ;  >  > 
else  < i ls*-dvl ; i2s=0; iend=60; > 
if  (dul>0)  {jlss0; j2s=dul ; jend=jrng-dul ; >  ng;> 
else  <jls“-dul ; j2s*0; jend=jrng; > 
if ( (vof«dvolay (vimg, ils,  i2s, iend,  j Is, j2s,  jend) ) (minvofs) 
<minvofs“vof ;mdu=dul ;mdv*dvl ; > 
printf ("t“*3d  n**3d  dul^Sd  dvl=X3d  vof=%lf  \n", 
t, n, dul, dvl, vof ) ;  >  /*n*/ 
for  <n*l jn (*t ;n++)  <  dvl — ; 

if  (dvl>0)  ■£  i  ls=0;  i2s=dvl ;  iend=£0-dvl ;  >  > 
else  < i ls“-dvl ; i2s=0; iend=60; > 
if  (dul)0)  { j ls=0; j2s=dul ; jend-jrng-dul ; >  ng ; > 
else  -C  j  ls*-dul ;  j2s=0;  jend=jrng ;  > 
i f ( (vof=dvolay (vimg, i 13, i2s, iend, J Is, j2s, jend) ) (minvofs) 
■Cminvofs^vof  ;mdu=dul  ;mdv*=dvl  ;  > 
pr i nt f ( " t  =X3d  n**3d  dul=%3d  dvl=%3d  vof*%lf\n  ", 
t, n, dul, dvl, vof ) ; >  /*  n  */ 
if  (mdv)dvl  +  l  &&  mdu(dul-l)  f  1  =  1 ;  else  f2>=0; 
printf  ( "mdu*X3d  mdv»*3d  fl»%2d  f2*%2d\n",  mdu,  mdv,  f  1,  f2)  ; 
if  (fi—i  f2«l)  break; 
for  (n*l ;n (*t+l ;n++)  <  dul — ; 

if  (dvl)0)  <  i  ls**0;  i2s*»dvl ;  iend«S0-dvl ;  >  > 
else  < i ls*-dvl ; i2s«0; iend=60; > 
if  (dul)0)  <j ls»0; j2s=dul ; jend”jrng-dul ; >  ng ; > 
else  <jls*-dul ;j2s*0; jend=jrng; > 
i f ( (vof»dvolay (vimg,  ils,  i2s,  i en*,  j Is,  j2s,  jend) ) (minvofs) 
<minvofs»vof ;mdu=dul ;mdv=dvl ; > 
pr i nt f  ( " t *X3d  n=*3d  dul**%3d  dvl=X3d  vof**lf\n  ", 
t, n, dul, dvl,  vof ) ; 

>  /*  n  */ 

for  (n=l ;n <=t  +  l ;n++)  <  dvl++; 

if  (dvl)0)  < i ls=0: i2s=dvl ; iend=60-dvl ;  >  > 
else  { i ls*-dvl ; i2s=0; iend=60; > 
if  (dul>0)  <jls»0;  j2s=dul ;  jend=»jrng-dul ;  >  ng;> 
else  -Cjls—dul ;  j2s=0;  jend-jrng;  > 
if ( (vof=dvolay (vimg, ils,  i2s,  iend,  jls,  j2s,  jend) )  (minvofs) 
■Cminvof s»vof  ;mdu=dul  ;mdv=dvl ;  > 
printf ("t=X3d  n«X3d  dul^fcSd  dvl=%3d  vof=%lf\n  ", 
t, n, dul,  dvl,  vof) ; 

>  /*n*/ 

if  (mdv(dvl-l  &&  mdu)dul+l)  f2=l;  else  fl=0; 

printf  (  "mdu=*3d  mdv=*3d  fl=*2d  f2=*2d\n" ,  mdu,  mdv,  f  1 ,  f2)  ; 
if  ( f  1  ==1  &&  f2=>=l)  break; 

>  /*t*/ 

d u C 1 3 =mdu ; dv C 1 1 =mdv ; 

>  /*  layer  */ 
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makevirng (ofs0,  ofsl, 1, zimg0, zimgl,  vimg,  cover) ; 

if  <ofs0)ofsl)  jrng«150-ofs0;  else  jrng*150-ofsl ; 
minvofs*dvolay (vimg, 0, 0, 80, 0,0,  jrng) ; 

printf  ("iiiinvofs*Xlf  dul*X3d  dvl*X3d  jrng*X3d  cover*Xlf\nH, 
minvofs,  dul,  dvl,  jrng,  coverCll ) ; 
for  (t*l ;t (30;t+*2)  < 

for  (n*l ;n (»t ;n++)  <  dul++; 

if  (dvl>0)  < i ls=0 ; i£s=dvl ; iend=60-dvl ;  >  > 
else  < i ls*-dvl ; i2s=0; iend=60; > 
if  (dul>0)  < j ls*0; j2s*dul ; jend^jrng-dul ; >  ng;> 
else  <jls«-dul ; j2s*0; jend“jrng ; } 
if <  (vof«dvolay (vimg, i Is, i2s,  iend,  j Is, j2s, jend) ) (minvofs) 
<minvofs«vof  ;mdu=»dul  ;mdv*dvl ;  > 
printf (Mt-X3d  n-*3d  dul«X3d  dvl«X3d  vof*Xlf  \n", 

.  t, n, dul, dvl, vof ) ;  >  /*n*/ 
for  (n*l ;n(«t ;n++)  <  dvl — ; 

if  (dvl>0)  ■£  i  ls=*0;  i2s=dvl ;  iend»60-dvl ;  >  > 
else  <ils»-dvl ; i2s=0; iend =80; > 
if  (dul>0)  <jls«0; j2s*dul ; jend=jrng-dul ; >  ng ; > 
else  < j ls*-dul ; j2s*0; jend=jrng ; > 
if ( (vof»dvolay (vimg, ils, i2s,  iend,  jls,  j2s,  jend) ) <minvofs> 
<minvofs*vof ;mdu*dul ;mdv=dvl 5 > 
printf ("t»X3d  n*X3d  dui«X3d  dvl*X3d  vof*Xlf\n  ", 
t, n, dul, dvl, vof > ; >  /*  n  */ 
if  (mdv>dvl  +  l  44  mdu(dul-l)  fl=»l;  else  f2*0; 
printf ("mdu»X3d  mdv«X3d  fl*X2d  f2*X2d\n", mdu, mdv, f 1, f2> ; 
if  (fl««l  44  f2**l)  break; 
for  (n*l  ;n (*t  +  l  ;n++)  <  dul  — ;. 

if  (dvl>0)  < i ls*0; i2s*dvl ; iend*60-dvl ;  >  > 
else  <i ls*-dvl ; i2s*0; iend«60; } 
if  (dul>0)  {j ls*0; j2s*dul ; jend*jrng-dul ; >  ng;> 
else  <jls»-dul ; j2s»0; jend*jrng; > 
if ( (vof«dvol ay (vimg, ils, i2s,  iend,  j Is,  j2s, jend) ) (minvofs) 
•Cminvofs*vof  ;mriu*dul  ;mdv*dvl ;  } 
printf ("t»X3d  n*X3d  dul*X3d  dvl=X3d  vof*Xlf\n  ", 
t, n, dul, dvl, vof) ; 

>  /*  n  */ 

for  (n*l ;n (=t+l ;n++)  {  dvl++; 

if  (dvl)0)  < i ls*0; i2s*dvl ; iend=60-dvl ;  >  > 
else  < i ls*-dvl 1 i2s"0; iend=60; > 
if  (dul)0)  < j ls®0; j2s»dul ;  jend=*jrng-dul ; >  ng;> 
else  <j ls»-dul ; j2s*0; jend=jrng ; > 
if ( (vof=dvo lay (vimg, ils, i2s,  iend,  j Is,  j£s,  jend) ) (minvofs) 
•Cminvofs*vof  ;mdu=dul  ;mdv=dvl ;> 
printf ( "t“X3d  n-X3d  dul*X3d  dvl»X3d  vof*Xlf\n  ", 
t,  n, dul, dvl, vof) ; 

>  /*n*/ 

if  (mdv(dvl-l  44  mdu)dul+l)  f2»l;  else  fl=0; 

printf ("mdu»X3d  mdv=X3d  fl=X2d  f2=X2d\n", mdu, mdv, f 1, f2) ; 
if  (fl**l  44  f2*“l)  break; 

>  /*t*/ 

d u C 1 1 =md u ; d v C 1 3 *md v ; 

>  /*  layer  */ 
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blcnt=fread ( (char  *) zofs0, 12C0, 1, in) ; 
fclose(in) ; > 

else  print f < "fi leopen  rcp\zofs03  failed\n"); 
return; > 

zofslinO 

</*  recovers  offsets  stored  by  zofsloutOin  file  *rcp\zofsl3*  */ 
int  blent; 

FILE  # in ; 

i f ( ( in=fopen ( "rcp\zofsl3", "rb" ) ) ! =NULL>  < 
blent=fread ( (char  *) zofsl, 1200, 1, in) ; 
fclose (in) ; > 

else  printf ( "f i leopen  rcp\zofsl3  failed\n">; 
return; > 

print  zofs ( ) 

{  int  i ; 

printf ("zofs  t=0\n"); 
for  (  i=0;  i  (150;  i++)  printf("X3d  %91f 
printf < "zofs  t=l\n"); 
for  ( i  =0 ;  i  <  150 ;  l  ++ )  printf("*3d  X91f 
ret  urn; > 

screenzofsO 
C  int  i,j=0; 
printf ("in  screenzofss key\n" ) ;getchar() ; 

INI TBC  ( )  ; SCREEN  (1,  1,0,0)  ; 

for  (i«0;i  (150;i++)  PSET (2*i, 150- ( int > < zofs0Ci 3 > , 1 ) ; 
for  (i=0;i (150;i++)  PSET (2*i, 200- ( int >( zofsl C i3 ), 2) ; 
printf ("*") ; get char ( > ; SCREEN (2, 0, 0, 0) ; SCREEN (0, 0, 0, 0) ; 

> 

printdwsO  < 
int  i ; 

for  ( i=*0;  i  <nl ;  i++) 

printf  ( "5(3d  dw0  *3d  dwl  X3d  zdw0  %3d  zdwl  X3d  du  %3d  dv  X3d  \n", 
i,  dwCil C03,  dwCiD [13, zdw[i3 C03 , zdwC i 3 C13 , duCi3 , dvCi3 > ;  > 

printvirng  ( jrng)  int  jrng;  <  /*  int  vimg Ct 3  C i 3  Cj  3  */ 
int  i , j , t 1 , t2; 

printf("\n  vimg  row  totals  for  jrng=*X3d\n" ,  jrng )  ; 
for  ( i=0; i (60; i++)  < 

for  (j=»0,  tl*0,  t2*0;  j  <10;  j++) 

■Ct  l+=vimg  [03  C  i  3  C j 3  ; t2+=vimg C 1 3  C i 3  C j 3  ;  > 
printf ("\n*6d,tl,X6d,t2\n",tl, t2) ;>  > 

/*  **************************************************************  */ 

Compiled  functions 
FEBLIB.C 

(♦include  "stdio.h" 

♦♦include  "math.h" 


",  i,  zofs0Ci3 ) ; 
",  i,  zofsl Ci3 )  ; 
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void  grab<a, b, c)  unsigned  a, b, c;  < 
out p (0x301, 0xac) ; 

while  (( inp (0*300)  ft  2)«=2)  continue; 

while  <(inp(0x300)  ft  2)**0)  continue; 

outp(0x30f,  a)  ;  outp(0x30e, b) ;  outp(0x30d, c) ; 
while  ((inp(0x300)  ft  2)»2)  continue; 

while  <(inp (0x300)  ft  2)*=0)  continue; 

out p (0x301, 0x2c) ; 
return;  > 

void  copyzimg  ( img0,  zimg0,  zimgl,  t,  q,  b) 

/•copy  img0  to  zimg I t=0, 1 lq*[23  i“[603  j»C1503*/ 
char  *irng0; 

int  zimg0C23 [60] [1503, zimgl [23 [603 [1503 ; 
int  t,q,b; 

<  int  i, int, j, jm, a; 

for  < i  =  <q-l ) *30, im=0; i <q*30; i++,  im+*720)  < 

for  < J=0, jm*30+im; j <150; j++, jm+=2)  <  a»  <  int)  img0[ jm3  ft  255; 
if (t*=0)  zimg0Cb3 Ci3 Cj3*a;  else  z imgl Cb3 Ci3 [j3*a; 

>  /*j*/ 

>  /*i*/ 
ret  urn ;> 

void  zolay0(zimg0, zofs0)  int  zimg0C23 C603 C1503 ;  double  zofs0C1533; 

<  int  ofs, i, j ; 

double  sumsq, meena»0, meenb“0, dmeen, del ; 
for  ( i=0; i (60; i++)  < 
for  j=«0;  j  (150;  j++)  < 

meena+* (double) zimg0[03  Ci3  C j 3 ;meenb+“ (double) zimg0C13  C i 3 Cj3 ; > 

}  dmeen=meena/meenb; 

printf ( "meena, meenb, dmeen  %S.  If  %9.lf  *9. 21 f\n", meena, meenb, dmeen) ; 
for  ( i*0; i (60; i++)  <  for  ( j«0; j (150; j++)  zimg0C13 Ci3 [j3**dmeen; > 
for  (ofs*»0;ofs <150;ofs++) 

<  for  ( j=0, *umsq*0; j ( ( 150-of s) ; j++) 

•C  for  ( i*0;  i  (60;  i++) 

•C  del=zimg0C03  Ci3  [j3-zimg0C13  Ci3  [j+ofs3  ;  sumsq+»del*del ;  >  > 
zofs0[ofs3*sqrt (sumsq/ ( 150-ofs) ) ; 

printf ( "zofs0CX3d3=*91f  sumsq*?^. 1 f  \n  ", ofs, zofs0Cofs3 , sumsq) ; 

> 

return;  >  1  . 

void  zolayl (zimgl, zofsl )  int  z img 1 [23 C603 C 1503 ;  double  zofsl[1503; 

<  int  ofs, i, j ; 

double  sumsq,  meena=0, meenb“0, dmeen, del ; 
for  ( i=0; i <60; i++)  < 
for  ( j«0; j <150; j++)  < 

meena+«  (double)  zimgl  [03  [i3  [j3  ;meenb+=»  (double)zimgl[l][i3[j3;> 

>  dmeen=meena/meenb; 

printf ( "meena, meenb,  dmeen  *9. If  *9. If  %9.  21 f\n", meena, meenb, dmeen) ; 
for  ( i=0;  i  <60;  i++)  <  for  (j=0;  j  <150;  j++)  zimgl  [13  [i3  [j3*=»dmeen;  > 
for  (ofs=0;ofs <150;ofs++) 

<  for  ( j=0,  sumsq=0; j < ( 150-ofs) ; j++) 

{  for  ( i=0; i (60; i++) 

{  del=zimgl [03 [i3 [j3-zimgl [13 [i3 [j+ofs3 ;  sumsq+*del*del ; >  > 
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zofsl Cofs3«sqrt (sumsq/ ( 150-ofs) ) ; 

printf ("zofsl C*3d3-X91f  sumsq-%9.  21  f  \n  ”,ofs,  zofsl  Cofs3,  sumsq)  ; 

>  return;  > 

void  makevimg  <ofs0,  ofsl,  1,  zimg0,  zimgl, vimg,  cover) 
int  ofs0, ofsl,l; 

int  zimg0C23  C603  C1503,  zimgl  C23  C603  C1503, vimgC23  C603  C1503  ? 
double  coverC103; 

C  int  i=0, J=®» zerolev,  tml,  tm2; 

double  sigma, sumsq =0, del, den, meana^0, meanb-0, dmean; 
for  <  i=0;  i  <60;  i++)  ■£ 

for  (j-0; j <150-ofs0? j++> 

<meana+-zimg0C03 Ci3 Cj3 ;meanb+«Timg0C13 Ci3 Cj+ofs03 ; >  > 
dmaan-meana/meanb ; 

printf ("X9.21f  %9.21f  *9. 21f \n", mcana, meanb, dmean) | 
den-1 . 0/ (60* ( 150-ofs0> ) ; 
zerolev*  meana*den; 

printf ( Mden=*lf  zero-level-%5d\n",  den,  zerolev) ; 

/*  calculates  arid  sets  zero  level  for  points  not  in  this  layer.*/ 
for  < i*0; i <60; i++)  -C 

for  ( j=150-ofs0; j <150; j++)  vimgC03 Ci3 Cj3«zerolev; 
for  < J-0; J <150-ofs0; j++)  < 

del-zimg0C03  Ci3Cj3-<int) (zimg0Cl 3 C i 3 C j+ofs03*dmean)  ; 
sumsq+-del*del ; >  > 
sigma-sqrt  <sumsq*den) ; 

printf  ( "ofs-X3d  sumsq-X9.  01  f,  siyma-*9.  21f\n", ofs0,  sumsq,  sigma)  ; 
for  < i*0; i <60; i++)  < 

for  < j*0| j <150-ofs0| j++)  < 

del«zimg0C03  Ci3  Cj3-(int>  (zimg0C13  Ci3  Cj+ofs03 "dmean)  f 
if  ( (del-fabs (del) ) (sigma)  <vimgC03 Ci3 Cj3»zimg0C03 Ci3 Cj3 ; 
cover C13++I J 

else  vimgC03 Ci3 Cj3-zerolev; 

>  > 

cover  C 1 3 "cover  Cl  3  *den ; 
sumsq-0;meana-0; meanb-0; 
for  (i*0; i <60; i++)  < 
for  (j"0; j <150-ofsl ; j++) 

•Cmeana+«zimgl  C03  Ci3  Cj3  ;meanb+»zimgl  C13  Ci3  Cj+ofsl3  ;  >  > 
dmean-meana/meanb ; 

printf <"X9. 21 f  %9.21f  X9. 21 f\n", meana, meanb, dmean)  ; 
den-1 . 0/ <60* < 150-of si ) ) ; 

zerolev«meana*den;  print f < "den-%1 f  zero-level-X5d\n",  den,  zerolev) ; 

/♦calculates  and  sets  zero  level  for  points  not  in  this  layer.*/ 
for  < i=0; i <60; i++)  < 

for  < j-150-ofsl ; j <150; j++)  vimg C 13 C i3 Cj3 -zerolev; 
for  <j-0; j  < 150-of si ;  j++>  < 

del -zimgl C03  Ci3  Cj3- (int) < zimgl Cl  3  Ci3 Cj+ofsl3*dmean) ; 
sumsq+-del*del ; >  > 
sigma-sqrt <sumsq*den) ; 

print  f  (  "ofs*X3d  sumsq»X9.  01  f,  si gma-*9.  21  f \n",  of  si,  sumsq,  sigma)  ; 
for  < i=0; i <60; i++)  { 
for  ( j-0;  j  <150-ofsl ;  j++)  ■( 

t ml -zimgl C03  Ci3  Cj3 ;t m2- zimgl C 1 3  C i 3  Cj+of si 3 ; 
del-tml-(int) <tm2*dmean) ; 
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/♦printf ("ij*%3d  %3d  zimgl  [03*X5d  zimgl  [13=X5d“,  i,  j,  tml,  tm2)  ; 
printf ("tmSadmean^^lf  (int)<  .  )*X5d,  del=*Xlf", 

tm2*dmean,  (int)  (tm2*dmean) ,  del )  ;  */ 
if  ( <del*fabs (del ) ) (sigma) 

vimgCID C i 3 Cj3*zimgl C03 C i 3 Cj3 ; 
else  vimgtlD Ci3 Cj3*zerolev; 

/•print f < "%4. 21 f  %3d  %3d  ",  del , z img 1 C03 C i 3 C j3 , vimg C 1 3 Ci 3 C j3 ) ; */ 
>  >  ret  urn ;  > 

double  dvolay (vimg,  i Is,  i2s,  lend,  j is, j2s, jend) 
int  vimg C23 C60D C 1503 ; 
int  i Is, i2s, iend, jls, j2s, jend; 

<  double  smsq*0; 

int  il, i2, jl, j2,  k*0,del ; 
for  ( i l*»i  Is,  i2*i2s;  i l  (iend;  i  l+-» ,  i2++)  < 
for  (jl=*jls,  j2»j2s;jl  (jend;jl++,  j2++)  < 

del=vimg C03 Ci ID Cj 13-vimg C 1 3 Ci2D Cj23 ;  smsq+=del*del ;k++;  > 

> 

smsq=smsq/k;  ret  urn  <  sr.sq )  ; 


void  showclouds(img0, zimg0, zimgl) 
char  *img0; 

int  zimg0C23 C603  £1503, zimgl C23  £603  £1503  ; 

< 

/*  converts  zimgt Eb3 Er3 Ec3  to  chars  and  writes  to  disk#/ 

/«  then  waits  for  svpcip  command  to  display  the  file  #/ 
int  q,  i,  im,  j,  k, blent ; 

FILE  *out,  *fopen(); 

if  ( (out»fopen( "rcp\zimg3. 000",  Hwb") )**NULL)  < 

printf <"f ileopen  rcp\zimg3.  000  fai led s key\n" ) ;getchar() {exit () {> 
for  (k-0;k <21600; k++>  img0EkD«0; 
for  <q=*0;q <2;q++)  < 

for  ( i»q*30,  im*i*360; i < <q+l ) *30; i++,  im+»*360) 

</*lst  q,  both  b’s,  t«l*/ 
for  ( j*0, k*im; j <150; j++,  k++)  < 
img0Ck+303  =  (char) zimg0E0D  Ei3  EjD ; 
img0Ck+ 1803* (char) zimg0£l 3  Ci3Cj3;>>> 
blcnt®fwrite( (char  *) img0, 21600,  1,  out ) {  /*  ql,2;  tl;  bl,2  */ 
for  (q=0;q <2;q++)  < 

for  ( i=q*30,  im«i*360;  i  <  (q+1 )  *30;  i++,  im+**360)  < 
for  ( j*0, k*im; j <150; j++,  k++)  < 
img0Ek+30D* (char) zimgl£0D£i3Ej3; 
img0Ck+ 1803=*  (char)  zimglE13CiJEjD;>>> 
blcnt=fwrite ( (char  *) img0, 21600, 1, out ) ;  /*  ql,2;  t2;  bl,2  */ 
fclose (out ) ; 
ret  urn ; > 


LINKER  TO  FEBLIB.C 


This  1 ibrary, FEBLIB,  contains  functions  for  WHSKY. 200 
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30“program  size 
2*memory  needed 

void  grab (unsigned, unsigned, unsigned) 

/•sets  address  for  grabber  ram  */ 
void  copyzimg (*,*,*, int,  int,  int ) 

/•copies  img0  into  zimg0, 1  to  be  filed*/ 
void  zolay0(*,*> 

/•days  zimg0C0ijl  onto  ,zimg0Cl  i  j3,  returns  zofs0*/ 
void  zolayl (*,  *) 

/•olays  zirngl  C0i  j"<  onto  zimglClijJ,  returns  zofsl*/ 
void  makevimg ( int, int, int, *,*,*, *) 

/•copies  1  lr  of  zimg0, 1  into  vimgCtl*/ 
double  dvolay (*, int, int, int, int, int, int) 

/•olays  vimgCts03  onto  vimgCt»13*/ 
void  showclouds(*, *, *) 

/•makes  a  f i le, rcpzimg3. 000,  to  put  on  monitor*/ 


Annex  A4. 2  MYLIB.C 

# i nc 1 ude  " st  d i o. h " 

♦♦include  "math.h" 

void  fi Iter (xa, xla, nf ft , e, r, vin)  /*  by  Aubanel  i  Byte  2.84  */ 
double  *xa, *xla;  int  nfft,ej  double  *r, *vin;  < 

/*  enter  with  xa, ya  filled  from  imgx  with  size  and  no.  of  freqs  e*/ 
/*  fill  fft  series  into  r  and  vin  */ 
int  j,jl,  j2,  jpl,  jp,  n-nfft,  k,  kl,  lj 
int  q=0, d«n/10, esq“e*e$ 
double  pi*3.  141592654,  twopi-2*pi ; 
double  twopiovn»twopi/n,  sl=0,  s2>0,  xl,  x2,  dm,  b; 
double  g«0,  s,  c,  ul,  f 1, f2, t, lg2, f ; 
double  uC923,  vC923 ; 
lg2*log (2. 0) | 

for  <j“0; j <d$ j++)  <sl+*xaCj3;  s2+-xaCn-j-13 j > 

/*  map  to  straight  line  */ 
xl*sl/d{ x2**s2/d5dm* <x2-x  1 )  /  (n-d)  ;  b*  < xl+x2-dm/n)  /2j 
for  (ja=0;j  <njj++)  <xaCj3+*-dm*j-b;g+*xaCj3  ;  > 
r  C03*g/n;q=*l ; 

j2*(n-l)/2;  jpl»log ( (double) <2*j2-l )) /lg2;  /*r(k)  transform  */ 
for  (k*qjk  (e;k++)  <  jl=j2;s=twopiovn*k;ca>cos(s)  i5*sin(s)  ? 

for  (j*l;j(*jl;j++)  < l=2*j-l |uCj3“xaC13*c+xaCl+13 ;vCj3=xaCl 3*s; > 
s«=2*s*c ;  c»2*c*c- 1 ; 

for  ( jp=*l ;  jp  (“jpl ;  jp++>  <uCjl  +  13=0;vCjl  +  13»0;jl“(jl+l) /2? 
for  ( j*l ; j <=j 1 ; j++>  <l=2*j-l5  ul=uC13*c-vC11*s+uCl+13 ; 

vCj3ssuC13*s+vC13*c+vCl  +  13  ;  uCj3=ul|>  /*j*/ 
ss»2*s*cjc®2*c*c~l ;  >  /*jp*/ 
r Ck3"(xaC03+uC13*c+vC13*s) /nj  >  /*k*/ 
for  (k^qik (e;k++)  <  /*  vin(k)  transform  */ 
jl*j2;s=twopiovn*k5c=cos(s) ;s=sin(s) j 

for  ( J®1 ; j (=j 1 ; j++)  <  l=2*j-l ; uC j3  =-xa  C  13*s;  vtj3*xa  Cl  3  *c+xa tl+13  ;  > 
s»2*s*c ; c=2*c*c- 1 ; 

for  (jp=l ; jp<*jpl 5 jp++)  <uC j 1+1 3 =0; vEj 1+1 3=0; j 1= ( j 1+1 ) /2; 
for  (j*l J J  <=jl ;  J++)  < l=2*j-l ;ul=uC13*c-vC13*s+uCl  +  13 ? 


& 


Micro  Science  Inc 


Page  46 


AFOSR  F49620-86-C-0093 


vCj3»uC13*s+vC13#c+vCl+13 ;uCj3*ul ;  >  /*j*/ 
s*2*s*c  5  c«2*c*c- 1 ; >  /* j  p*/ 
vinCk3*-(uC13*c+vC13*s) /n;  >  /*k*/ 
if  (e>q)  q=e;  /*  inverse  transform  */ 

for  <k*l, f 1*0, f2*0;k<e;k++)  <  t*rCk3{  fl+*tj  f2+*k*k*tj> 
h laC03=r C03+2* (f 1-f 2/esq) +b} 
jpl*log ( (double) <2*e-3) ) /lg2; 
for  ( j  =  l ; j  <n; j++)  < 

for  (k=l ; k <e; k++)  {  f=l-k*k/esq;  uCk3=rtk3*f;  vCk3*-vinCk3*f ; > 
kl*e-l ;s*twopiovn*j ;c*cos(s) ;s=sin(s) ; 

for  (jp=l j jp<=jpl ; jp++)  {uCkl+13*0;vCkl+13*0;kl*(kl+l)/2; 
for  (k*l ;k <*kl ; k++)  <l*2*k-l ;ul*uC13*c-vtl3*s+uCl+13 ; 

vCk3*uC13*s+vC13*c+vCl+13 ;uCk3*ul ;  > 
s*2*s*c}c«2#c*c-l; >  /*jp*/ 
xlaCj3*rC03+2*(uC13*c+vC13*s)+dm*j+b|  >  /*j*/ 
ret  urn ;  > 


U  GENERAL  INFORMATION 


The  TN2505c»nera  series  uses  LSWLSJ  toad  state  construction.  enabling  the  camara  to  be  fight  hseigh*  and  compact  in  size.  These 

Advswceainthto  technology  have  produced  a  camera  that  has  low  cameras  are  ideal  lor  security.  ntBtary.  automation,  medical  and 

power  dissipation,  *id  whoa#  pertormancawB  not  degrade  over  sciertific  applications, 

time.  Solid  aUWdrcultryhaaaMrninatedbufcy  camara  components 


BOARD  CONFIGURATION 


A.  IMAGER  BOARD  —  Contains  the  CIO 
(Charge  Injection  Device).  LSI  scan  generator 
and  a  hybrid  video  pre-amplifier. 

B.  SCAN  GENERATOR  —  Located  on  tha 
lmage<'  board,  this  CMOS  scan  generator 
provides  the  dock  signals  required  by  the  CIO 
imager,  analog  video  processing  circuits  and 
tha  sync  and  blanking  pulses. 

C.  VIDEO  BOARD  —  The  video  amplifier  on 
this  board  receive*  video  from  the  hybrid 
pre-amplifier  (of  the  imager  board)  and 
provides  the  following  circuit  functions: 
automatic  black  level,  gain,  and  the  addition 
of  sync  and  blanking. 


D.  ROWER  SUPPLY  BOARD  -  This  board 
contains  the  DC  power  circuit  and  modular 
sync  component  system.  Power  is  supplied  to 
the  camera  by  a  separate  power  supply 
module.  A  lO*  power  supply  cable  is 
standard,  although  up  to  SO*  of  cable  may  be 
used  without  affecting  quality.  This  camera 
can  dso  be  operated  from  other  OC  power 
sources**)  output  of  +12  to  +33vDC  O 
1/2  amp  peak  minimum. 


PERFORMANCE  FEATURES 


OPTICAL  VERSATILITY — A 

Standard  C  mount  permits  the  use 
of  a  wide  variety  of  CC7V  lenses, 
InckxSng  telephoto,  zoom,  wide 
angle  dose-up  and  pm  hole.  The 
exact  application  would  determine 
which  lens  to  use.  An  Automatic 
Lens  Connector  provides  video 
and  power  tor  an  auto  iris  lens 
system. 


.r  SENSOR  QUALITY  —CIO 
seneon  exhibit  virtually  no  lag  or 
blooming.  Individual  pixel  image* 
•xn©n  nwwnai  OKXXTwig  into 
adjacent  pixels.  Thereto  no  poet 
Held  lag,  or  bufid  up  lag,  and  bias 
ttoht  rnmnuMtion  Is  not 
required. 


I/O  CONTROL 


OAMOPTKJNS  — Internal  vw*o  gain  of  INJECT  INHIBIT  —  Inject  Inhibit  allows  UMNO  OUTPUT*— Camera  timing 

+6dbrv +i2dbirtay  be  implemented  via  the  interruption  of  camera  readout  upon  signals  are  available  at  the  I/O  connector  and 

I/O  control.  When  the  camera  to  not  being  command,  synchronous  stop  motion  of  high  may  be  used  to  synchronize  the  operation  of 

used  with  an  auto  iris  lens,  an  optional  speed  events,  and  integration  of  static  low  the  extern*  equipment.  The  sensor  element 

external  gam  switch  can  be  adoed.  light  imagery.  When  Inject  Inhibit  to  removed,  rate  dock  sign*  is  provided  along  with  H&V 

a  single  few  of  information  may  be  stored  drive  signs*.  However,  these  signals  are  not 

externally  for  further  processing  or  monitor  buffered  and  require  care  in  Interfacing, 
display. 


SENSOR  ARCHITECTURE 


TheTN  2505  series  imagers  conatt  ola 
matrix  of  photo  sensitive  pbcete  arranged  In 

(columns  and  row*.  Tha  readout  it  condoled 
by  ahift  registers  (or  decoders)  located  at  the 
periphery  of  tie  imager.  Photon  reSation 
imaged  on  the  sensor  surface  penetrates  the 

I  sensor  surface  where  it  ia  collected  by  tha 

charge  storage  areas  of  each  pixat.  (Pig  1) 

Tha  column  storage  area  of  each  pixel,  aa 
we«  as  the  interconnection  of  each  column's 

I  storage  areas  for  pixels  in  that  column  ia 

ordered  by  a  single  strip  of  high  conductivity 
pofysWoon,  which  vertically  croasee  tha  thin 
oxide  regions  of  a>  pixels  in  a  column  of  the 

1CX)  matrix.  The  row  storage  area  iasimtariy 
ordered  by  a  horizontal  crossing  of  the  CIO 
matrix.  The  column  and  row  storage  areaa 
are  dWectrically  separated. 

8  A  thin  aluminum  strip,  in  c%ect  contact  with 

the  upper  potysMoon,  is  designed  to  improva 
imager  noise  performance  imd  row  access 
speed. 

The  charge  coflected  ijndsr  a  spedfk.' 

S  pixel's  electrode  area  can  be  transferred  or 
infected  between  column  and  row  of  (he  pixel 
storage  area  Thie  tonct  on  results  in  an 
imager  with  no  spacing  bsc.«een  pixels  and 

I  minimum  obscuration  to  Vicktant  Rumination. ' 
AddHionafty,  the  array  structure  to  “tuned" 


FIGURE  1 

to  minimize  multiple  reflection  Interference. 
Hence,  high  quantum  efficiencies  (45-55%) 
and  responsiveness  (approximately  0.2 
ampa/watt)  are  obtained, 

Rgur*  2  thowi  ttw  incorporation  of  t 
msntx  of  four  rows  and  four  columns  fid 
ptxsis  total).  These  form  a  differential  row 
read  CJO  imager.  Specific  columns  and  rowe 
are  addressed  through  column  and  row 
multiplexers  by  logic  r  output  from  the 


RGURE2 

column  and  row  selection  logic. 
SbnuRaneous  selection  of  a  column  and  tow 
causes  sigrtaf  charge  codec  ted  undor  tha 
cohetwi  MCifoos  vo  m  tiirwvrroQ  to  ro  row 
electrode.  There,  H  ie  sensed  and  ampUfled 
bye  high  gain,  low  notes  preamplifier.  Imager 
KTC  noise  la  effectively  sSminated  upon  tow 
selection  by  tha  currant  mode  amplifier. 
Column  induced  fbsed  panom  notes  ia 
removed  at  the  Imager  output. 


Tha  TN  25C5  camera  la  deal^ped  for  RS 170  output 
and  tha  TN  2506  for  COR  standard*. 


Pixel 


Active  Elements 


Size 


TN  2505  (CIO-1 


248V  x  388H 


242V  x  377H 


2 


TN  2506  (CIO-1 


294V  x  418H 


288V  x  400H 


TN2906A,  TN2506B 


SYSTEM  INTEGRATION 


Video  switching  or  imege  processing  requires  a  high  degree  of  sigrMl 
synchronization.  In  order  to  meet  these  requirements  the  TN2505  & 


A-1  UNX  LOCK  SYNCHRONIZATION  Une  lock  Sync  locks 
the  vert Icel  blanking  Interval  of  the  cemera  to  the  60  (or  50) 
cycle,  power  line.  When  several  cameras  ar*  used  Line  Lock  Syne 
eliminates  vertical  roll  when  switching  from  one  earner*  to  another. 
Unleee  specified  otherwise,  both  the  TN2S05 &  TN2006  come  with 
Une  Lock  Sync 

A-2  CffY*TAL  COMTROUXD  SYNCHRONIZATION  Crystal 
Controlled  Sync  provides  optimum  sync  stability  of  an  indhridual 
camera  operating  independently.  Tha  dock  crystals  supplied  tor 
these  cameras  are  TN2S05A  2-14. 318  MHz,  TN2S06A  2-15. 375 
MHz 

A-3  OSNLOCX  SYNCHRONtZATTOM  Genlock  Synchronization 
permit*  cameras  to  be  locked  to  each  other  or  locked  to  a 
separate  source  of  RS-170  or  CCtR  video.  These  models  have  an 
easy  access  Phase  Thm  Control  vdtoh  bringe  tha  Sync  Master 


TN2506  model  vsrdone,  A-1  through  A-d,  are  offered  with  a  variety 
of  eync  optlone  to  meet  virtually  any  situation. 

cwnera  and  tha  other  cameras  connected  to  the  master  Into  exact 
phase  relationship. 

A-4  H  &  V  DIWI  LOCK  SYNCMNOMZJKTION  This  version 
permits  one  or  more  cameras  to  be  locked  to  a  negative  going  H  d 
V  drive  of  an  RS-170  or  CCtR  video  source. 

A- SCOMRULOCK™  SYNCHRONIZATION  TheCompulock 
Sync  permits  one  or  more  cameras  to  be  locked  to  an  external 
source  of  progressively  teamed  (60/50  frames  per  eicond) 
video  auch  es  a  CRT  controller. 

A-d  COMRULOCX  ™ SYNCHRONIZATION  This  version  permits 
one  or  more  cameras  to  ba  locked  with  H  A  V  drive  outputs  of  a 
progressively  scanned  video  source.  Both  tha  A-5  and  tha  A-d 
versions  art  scanned  at  60  frames  par  second.  Whia  video  can  be 
dbpiayed  on  a  standard  monitor,  the  picture  w«  not  be  Interlaced. 

The  monitor  w*  lock  into  a  field  at  random. 
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U.  SPECIFICATIONS 


ELECTRICAL 

Input  VWMQt 
Power  Supply 


Camara 


Sanaor  Formal 


CompoaNe 

1' 


_ _ _  ALC  Range 

ttCMtC*  t0H<TN2SO5)  QrayScala 

50/BOHz  - 

220M»C±10H<TN2S06)50HZ  Modulation 

♦  ttjo+MVgC _  "franefer 

_ Fynctton — 

iLLltWtlr-at  SCOn  ROIO 


244  Varik*  306  Horizontal 
(TN2S09) 

nO  Vartical.  410  Horizontal 


1  V)tt  D-0  into  75  Own*  RS-170 


(Mama I  RS-170,  Unetock 


Sanaitkrtty 


300 10 1100  Nanometer* 


110  Com. 


2.0  MMon  to  1  WHhAO-300  EMVWONMEMlAL _ 

■»*>±!g.Lgl* -  Tamo.  Rang*  OjojOC. 

ttehadeaolgray.MMmtaii  HumkOhr  OSH  Non 

BA  Standard  Owl -  £3^ - 


Bal  tar  than  00Hal2S0anaa 
525/80  TN  2505 

62530  TN  2500 _ 

Ful  output  (1  Vjtt  ootnpoerte)  at 
,8<c  facanlata  (Ode  gain) 

Ate  faceplate  ( ♦  6 cm  gain) 
UaaaOia  vidoo  at  21c  facaplata 


♦  66  ♦  12dBGain, Sync, 
Infection  Inhfclt. 

♦  I2to  ♦  35VDC 


Video 

Cor  mac  Ion 


50Q(Yt  ekwwavaaMOme 

duration; _ 

OQ'a at 500Hz.  For appNcatbne 
requiring  graaiar  toadtog  w* 
recomma.iduaingtti*  two  piece 
camara  configuration. _ 


MECHANICAL 


iwatMwm 


For  further  Information  contact! 

Marital  Ing  Manager,  Bactronle  Camara  Operation 
GENERAL  ELECTRIC  COMPANY 
WO  7th  North  Street 
Liverpool,  New  Ytorit  134W 
(3  IS)  454-2SO 1/2945 

We  bring  good  things  to  life. 

GENERAL  HI  ELECTRIC 


*3  PWNTED  IN  USA 


fail*) 


I. 


Introduction 

The  TN2S0S,  6,  7  of fora  feature*  and  porfonaanco  aovor  boforo 
available  in  a  aurvolllanco  product.  Tha  camera'*  amall  alto  opona 
a  whole  realm  of  exciting  new  application  and  inatallation 
poaaibilitiea.  Ita  low  power  and  wide  input  voltage  range  further 
extend  ita  flexibility.  Uae  of  a  ac.  id  state  laager,  CE'a  CIO 
(Charge  Injection  Device),  eliainatea  the  performance  drift  end 
ultimate  need  for  replacement  experienced  with  vacuua  tube  imaging 
device* .  Utilixation  of  Large  Scale  Integrated  circuit*  (LSI' a) 
aakea  poacible  a  level  of  performance,  atability  and  reliability 
unattainable  with  dlacrete  component*. 

The  liat  below  outline*  the  many  feature*  which  make  the 
TN2S05,  6,  and  7  a  unique  aurveillanee  product. 


TM2303.  6.  7  Feature*  Liat 

e  Extremely  amall  alre  (3.37"x3.0"x3.04") . 

e  Light  weight  -  13. 5  ox.  Rugged  eonatructlon. 

e  Single  power  aource:  12-3S  VDC  at  2  watta,  or  117V  or  230 
VAC  50/60  cy  with  power  aupply  module. 

e  Standard  SS170  output  format  (TN2505,  7). 

e  Standard  CCIR  output  format  (TU2S06) . 

e  Available  with  line  lock,  or  cryatal,  or  ayne  lock  control. 

e  2.6  million  to  1  ALC  with  optional  automatic  irla  lena.  (fl.4 
to  f3(0) . 

e  Selectable  OdB,  6dB,  12dB  gain  booet. 

e  extremely  good  anti-blooming  performance. 

e  Inject  inhibit  feature  for  very  low  light  level  and 
"Snapahot"  application*. 

e  True  horizontal  aperture  correction  (enhanced  image 

aharpneaa) . 

e  Automatic  black  level  circuit. 

e  Adjuatable  vernier  mechanical  focua  Independent  of  lena. 

e  1/4x20  mounting  hole  on  top  and  bottom. 
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Tli*  following  specification*  represent  typical  camera 
performance.  General  Electric  Company  reserves  the  right  to  mat# 
changes  in  design  which  affect  specified  performance  without 
notification  or  obligation. 

IlfS.tgjgJl 

Input  Voltage  Range:  12.0  to  3S.0  OVDC 

Input  Power:  Slightly  over  2  watts  within  specified  Input 
voltage  range. 

Scan  Rates: 

SKSIAfa  LIN1L0CK 

Modal  Vertical  Horizontal  YacUtal 

TM2S0S , 7  S9.94HX  15.734HX  60Hn  15,730Hx 

TM2506  S0.00HZ  13.62SHX  SOHs  1S.62SHX 

Geometry:  Unaffected  by  input  voltage,  ambient  temperature, 
magnetic  fields. 

Signal-to-Noise  -  Peak  Signal  to  RMS  Noise:  SO  dB*  TM2S0S 

49  dB*  TM2506 
*at  Imager  Saturation. 

ALC  (with  Optional  Lens):  2.0  million  to  1 
Manual  Gain  Selection:  OdS,  <dS,  12dB 


Oisplayed  Resolution: 
(pixels) 

TM2S09 

Horlxontal  -  377 
Vertical  -  242 

CID17B 

TN2506 

Movlxontal  -  400 
Vertical  -  287. 5 

CXD18R 

TN2507 

Horlxontal  -  377 
Vertical  -  484 

CID20B 

Gray  Scale  Rendition: 

10  shades  of  gray,  minimum 
gray  scale  chart) . 

(El  A 

Sensitivity:  Pull  output  (IV  composite)  at  .8PC  faceplate 
illumination  with  OdB  boost  gain  selected. 

Pull  output  (IV  composite)  at  .027FC  faceplate 
illumination  with  injection  inhibit  period  of 
.5  second. 

Spectral  Response:  .35^  to  1.1  fJ  (see  Figure  1) 


Output  Format:  1.0V  composite,  Mbit*  poeltlva,  RS170 
(TU2505) ,  (TN2507) ,  CCIR  (TM2506) . 

Machanlcal 

Camara  Slza:  3.04"  long  (Laaa  Luna’ 

3.00"  wide 
3.37"  high 

Cmri  Walght:  13. 5  or. 


Camara  Mounting:  1/4-20  hola  top 

1/4-20  hole  bottom 


Mating  Connactora: 


(Viking  Induatrlea) 


Pwar  TKP-07-110-TS100T(Vikina  Induatriaa) 
Vidao  (BNC)  UG  1094/U 

Iria  TKP-07-110-TS100T  (Viking  Induatriaa) 


Temperature:  0°C  to  SO°C  I 

Relative  Humidity:  0-951  non  condensing 

Maxima  Altitude :  20,000  ft. 


[in?  nr? 


Electrical  (American)  (P/N  23289116) 

Input  Voltage:  120V  ♦  151 

$0/60  cycle,  2.S  watts 


Output  Voltage: 
Output  Current: 


24.0  OVpC 
180  HA  Mai 


noalnal 


Arc  Reference:  28  VRKS’ 


Maximum  Current:  1  MA 


Ilectrlcal  (European)  IP/M  232B9127] 

Input  Voltage:  220V  ±  101 

$0  cycle,  2.6  wa 


watts 


Output  Voltage:  24.0  OVDC  nominal 

Output  Current:  180  MA  jfai. 


AFC  Reference: 


28  VRHS 


Maximum  Current:  1  MA 


Mechanical  (American)  (P/N  232B9116) 

Size:  2.25"  long  x  2.25"  wide  x  1.62"  high 
Cable  Length:  124  ±  4" 

Weight:  11  oz. 


Size:  3.50"  long  x  2.12"  wide  x  1.65"  high 

! 

Cable  Length:  Input  36  ±  4"  Output  124  +  4 
Weight:  13  oz. 


Environmental 


Iifirttart:  0°C  to  S0°C 

Relative  Euaidlty:  951  bob  condensing 


■axlama  Altitude:  20,000ft. 


Any  duugi  incurred  to  equipeant  during  shipping  should  bo 
doeuaoatod.  It  Is  ths  responsibility  of  tho  carrier  to  insure  that 
equipment  arrives  is  serviceable  condition.  Mote  the  physical 
condition  of  the  packing  carton,  then  carefully  unpack  tho 
contents.  Inspect  each  item  for  daaage.  Damage  should  be  reported 
InsMdlately  to  the  responsible  agency.  Using  the  packing  slip  and 
the  list  of  equlpawnt  supplied  below,  ascertain  that  your  order  is 
correct.  Report  any  deficiencies  to  the  General  Electric  Coaqtany 
address  indicated  on  the  packing  list. 

IquipaMnt  Supplied:  TM2SOS,6,7  eaaera  without  lens 

Power  supply  Module  and  10  ft  interface 
cable. 

1/0  eonneetor  Viking  Industries  part 
nuaber  TKP-12-110-TS100T  with  12  pins 
(disassembled) . 


Power/Lens  Connector  (disassembled)  - 
Viking  Industries  part,  nuaber  TKP-07- 
110-TS100T  with  3  pins  (A2,A3,A4  swdels), 
5  pins  (A1  aodel). 
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SILICON  UIOEO  USER'S  MANUAL 

1.  INTRODUCTION 

Silicon  Uldeo  allows  a  PC  to  digitize,  procsss,  display,  and  transmit 
video  information.  Silicon  Uldso  will  digitize  ona  or  several  frames 
from  a  video  camera  or  other  video  source,  allow  the  PC  to  process  the 
image  data,  display  the  image  data  on  o  Btw  or  RSB  monitor,  and  allow 
the  PC  to  transmit  the  lease  data  via  a  modem  and  the  telephone 
network  to  a  similarly  equipped  PC  at  a  remote  location.  Uldeo  data 
is  digitized  and  displayed  at  8  bits  per  pixel. 

Silicon  Uldeo  can  contain  up  to  1MB  of  lease  memory  for  multiple  image 
storass  or  use  of  interlaced  imasss  greater  than  540  pixels  per  line. 

The  unique  features  that  dlstinsulsh  the  Silicon  Uldeo  from  other 
video  i easing  products  ore: 

Uarloble  sampling, 

1MB  image  memory, 

Uldeo  memory  address  rsgi stars, 

752  by  480  maximum  rssolutlon. 

Multiple  image  storage,  and 
Extensive  software  base. 

The  SUPCIP  software  (which  is  provided  ot  no  cost  for  each  Silicon 
Uldeo!  allows  the  user  to: 

Select  variable  resolution  sampling  parameters. 

Digitize  one  or  several  images. 

Display  one  or  several  images, 

Digitize  and  display  portions  of  the  raster, 

Rsod/wrlta  Images  or  portions  of  images  to/from  disk 
with  or  without  compression, 

Perform  addition,  subtraction,  differencing,  insertion, 
duplication,  and  averaging  on  two  Images  or  portions 
of  two  images. 

Perform  bit-strip,  histogram,  threshold,  exclusive -or, 
contrast  stretch,  reduction,  print,  half-tone,  and 
dither  on  ona  image  or  portion  of  an  image, 

Perform  convolution  or  dynamic  thresholding  using  a 
3  by  3  to  7  by  7  pixel  kernel  on  an  Image  or  portion 
of  an  image. 

Overlay  text  on  an  image  with  up  to  83  eharocters 
by  85  lines  available. 


Silicon  Uldeo  uses  conventional  dynamic  random  access  memory 
integrated  circuits  (DRAM 3 .  This  allows  memory  costs  to  be  kept  low 
and  allows  the  basic  design  to  increase  in  capacity  as  the  capacity  of 
DRAM  increases. 

Silicon  Uldeo  has  a  dual  ported  image  memory  that  allows  access  by  the 
PC  and  by  video.  While  the  image  memory  has  dual  access,  it  mow  not 
be  accessed  by  video  ond  PC  simultaneously.  The  SUPCIP  software 
accesses  the  image  memory  in  two  modes:  when  video  is  enabled  It 

accesses  during  the  vertical  interval,  and  when  video  is  disabled  it 
accesses  as  required.  This  allows  the  user  to  perform  image 
processing  and  view  the  results  as  the  processing  is  taking  place, 
while  providing  a  higher  speed  mode  when  repetitive  operations  ore  to 
be  performed  and  human  monitoring  of  the  processing  is  not  required. 


SILICON  VIDEO  USER'S  MANUAL 


1-3 

5*  sampling  and  display  of  pixels  is  control lad  by  a  micro  sequencer. 
This  allows  spscific  pixsls  to  bs  digitized  or  displaysd.  It  also 
sIIom  pixsls  to  ba  both  digitized  and  displaysd  on  ths  sons  roster, 
nts  programmable  portion  of  ths  micro  ssqusncsr  is  callsd  ths  XSAM  CX 
Sample  Memory 3 .  It  is  eonstuetsd  of  a  pixel  countsr  that  provides  ths 
•ddress  to  a  high  speed  memory  and  register.  At  the  start  of  each 
horizontal  line  ths  counter  is  cleared  to  zero.  At  the  trailing  edge 
of  horizontal  drive  (the  left  of  ths  raster!  ths  countsr  is 
incremented  by  each  pixsl  clock.  Ths  outputs  of  ths  memory  are 
clocked  into  a  register  that  controls  the  sampling  end  display  of 
Pixsls  and  other  functions.  Ths  XSAM  is  disabled  during  vertical 
blanking  or  when  the  PC  loads  the  bit  in  the  control  register  that 
disables  the  XSAM.  The  XSAM  addresses  that  correspond  to  ths  video 
roster  are  shown  below: 


Left  of  Raster 
I 

I  first  Nonblank  Pixel 

I  I 

i  I 

0  60  Active  Video 


Right  of  Roster 

I 

I 

I 

I 

BIS 


Note  that  there  is  a  maximum  of  7SS  pixels  that  are  active  video. 


The  XSAM  may  be  programmed  to  digitize  less  than  7SS  pixels  on  each 
ective  line  of  video.  The  video  may  bo  digitized  as  one  frame  of 
Information  (two  fields!  or  one  field  of  information.  Each  field  of 
Information  results  in  S40  linss  of  information.  Ths  highest 
resolution  la  obtained  using  a  frame  of  video  and  7SS  pixsls  per  lino. 
Nowsver,  this  rsquiros  the  1MB  image  memory  option.  The  maximum 
resolution  image  that  the  SS6KB  image  memory  will  support  for  a  frame 
ef  video  is  540  pixels  by  ths  460  active  lines.  Some  of  the  possible 
dumber  of  images  that  the  1MB  version  supports  erst 


6  Images  at  75S  pixels  by  480  lines, 

3  images  at  640  pixsls  by  4B0  lines, 

4  images  at  SIS  pixsls  by  4B0  lines, 

6  images  at  640  pixsls  by  S40  lines, 

B  images  at  SIS  pixsls  by  S40  lines,  or 

17  images  at  SS6  pixels  by  S40  linss. 


Nats  that  these  are  a  very  small  number  of  ths  possibilities. 


Your  PC*  Limited  28610  or  286*2  offer  *  number  of 
advanced  features  which  place  them  at  the  forefront  of 
powerful,  adaptable  yet  easy-to-use  personal  computers. 
Here  are  just  a  few  or  iuose  features: 

•  The  286s  work  at  both  the  original  IBM  PC-ATs  6 
MHz  meed,  and  also  at  either  10  MHz  or  12  MHz, 
depending  upon  which  model  you've  purchased. 

This  isn't  a  matter  of  just  changing  a  clock-crystal, 
but  rather  the  result  of  a  fresh  approach  to  computer 
design. 

This  increase  in  speed  translates  into  much  faster 
processing  of  data  in  your  work.  By  most  popular 
benchmarks  for  operating  speed  in  personal  com¬ 
puters,  the  28610  and  28612  are  up  to  ten  times  faster 
than  the  original  IBM  PC 

The  actual  increase  in  total  system  performance  you’ll 
find  when  working  at  your  286  (or  ‘‘throughput,’*  in 
computer  lingo)  will  depend,  on  the  programs  you 
use  and  the  speed  of  connected  peripherals,  such  as 
modems  and  printers. 

•  Your  PCs  Limited  286  is  dramatically  smaller  than 
earlier  designs.  By  saving  as  much  as  25%  of  the 
usual  “footprint,"  or  space  occupied  by  the  com¬ 
puter's  cabinet,  you  regain  space  on  your  desk  for 
ocher  items. 

•  The  six  16-bit  (PC-AT-style)  expansion  slots  and  two 
8-bit  (PC-style)  expansion  slots  in  the  PCs  Limited 
286s  give  you  more  for  your  money  by  allowing  you 
to  expand  your  computer  with  add-in  modems,  video 
cards,  multipie  I/O  ports,  expanded-memory  cards, 
speech-recognition  systems  and  many  other  exciting 
products.  Unlike  other  “AT  compatibles,”  the  PCs 
Limited  286s  don't  cripple  your  ability  to  expand 
your  computer  to  meet  your  changing  needs. 

•  As  you  know,  your  PCs  Limited  286  costs  less— 
often,  much  less — than  older,  slower  designs  with 
fewer  features. 
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Appendix  J 
Specifications 


au*,  large  chassis. . 

.  W*K3D  ‘*ji 63x53  3x4J$cm 

sae,  sn  ail  chassis ... . 


...  65  Hx  18.75  W  x  16 J  D  UU 165x47 6 x 41 9  cm 
. . SS  tbsf25  bias 

dock  speed  (user-selectable,  from  keyboard) . 

. . . 6  MHx,  10  MHx  (28610) 

8AM division  option  (user-selectable)... . . 

. * . Proirams>  silKfUr  VDISK 

. 840Kfor  DOS  programs,  384Kfor  VDISK 

. . . .  256K  for  DOS  programs 

. .  for  DOS  programs 

: : : :  PT‘"“ 

. . I  fidl-height,  3  hZf-height,  large  chassis 

. . 3  hay-height,  small  chassis 

expansion  bats  . . . 

. 8  IBM  PC-AT-compatible  (16-bit) 

"£***  (^sdeciabuZ. 

power  supply . . 

typical  power  consumption  at  llOv,  with . 

. .  V  •  “uu  'JJi  Mbflomhdisk  IS  warn 

. !•*  ™b,  360  Kb  floppy  disk  drives,  18  watts 

. ^  MbfiflPPy*  ^0 Mb  hard-diskdrive,  24  watts 

. 1-2  Mb,  40  Mb,  60  Mb  tape  drive,  33  watts 
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